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Recently,  a  description  of  the  molecular  mechanism  of  retention  in 
reversed  phase  liquid  chromatography  (RPLC)  was  derived  by  using  a 
partition  model  and  statistical  mechanical  theory.   This  theory  relates 
the  capacity  factor,  k',  of  a  solute  in  RPLC  to  the  binary  interaction 
constant  between  the  solute  molecules  and  the  mobile  phase  molecules. 
The  binary  interaction  constant  is  related  to  the  interaction  free 
energy  through  the  Boltzmann's  constant  and  the  absolute  temperature. 
Solute-solvent  interaction  free  energies  can  be  obtained  from 
chromatographic  data  if  the  partition  theory  of  solute  retention  is 
correct . 

We  have  tested  the  validity  of  the  predictions  from  this  theory  by 
using  more  than  300  sets  of  retention  data  from  the  literature  and  our 
laboratory.   With  only  a  few  exceptions,  the  data  confirm  all  'the 
predictions,  and  the  solute-solvent  interaction  free  energies  are  found 
to  be  solute  size  dependent,  which  was  a  case  not  considered  when  the 


theory  was  first  developed. 

Although  RPLC  is  one  of  the  most  popular  separation  techniques, 
the  actual  effect  that  the  stationary  phase  has  on  solute  retention  is 
still  not  completely  known.   Much  research  has  been  done  on  the  role  and 
effects  of  the  mobile  phase  on  solute  retention  in  RPLC,  but  there  has 
been  a  lack  of  systematic  approaches  in  characterizing  RPLC  columns. 

We  have  developed  a  simple  method  to  classify  RPLC  columns  in 
accordance  with  their  retentivity  by  using  more  than  20  solutes  of  very 
different  chemical  structures.    The  solvent  strength  of  the  mobile 
phase  has  to  be  accounted  for  before  the  retentivity  of  these  columns 
can  be  investigated.   This  is  accomplished  by  plotting  In  k '  of  the 

solute  against  the  ET(30)  polarity  values  of  the  mobile  phase.    A  good 
parameter  that  sums  all  the  retention  due  to  the  stationary  phase  is  the 
capacity  factor  of  the  solute  at  100%  water.  In  kw.   By  utilizing  In  kw 
and  the  slope  from  the  ET(30)  plot,  a  relative  retentivity  of  RPLC 
columns  is  obtained. 


CHAPTER  I 
INTRODUCTION 

Reversed  phase  liquid  chromatography  (RPLC)  is  the  most  common 
mode  of  high-performance  liquid  chromatography  used  today.   RPLC  is 
estimated  to  account  for  more  than  80%  of  chromatographic  systems 
currently  employed  (Melander  and  Horvath,  1980) .   The  name  "reversed 
phase"  was  introduced  by  Howard  and  Martin  (1950)  when  they  used 
non-polar  liquid  paraffin  and  n-octane  as  the  stationary  phase  to 
separate  fatty  acids.   They  used  the  term  to  distinguish  it  from  "normal 
phase"  chromatography,  which  consists  of  a  polar  stationary  phase  and  a 
non-polar  mobile  phase.   The  use  of  RPLC  has  shown  an  overwhelming 
increase  when  columns  packed  with  chemically  stable  microparticulate- 
bonded  silica  became  available.   The  practical  use  of  commercial  bonded 
phase  packings  in  RPLC  can  be  dated  back  to  Kirkland  and  DeStefano 
(1970) . 

These  bonded  phases  add  great  advantages  to  RPLC,  including  the 
fact  that  aqueous  mobile  phases  having  low  toxicity  and  high  optical 
transparency  can  be  adopted  to  achieve  most  separations.   Also,  the 
stability  of  these  bonded  phases  are  superior  to  other  stationary 
phases.   Aqueous  mobile  phases  of  pH  between  7.5  and  2.5  are  compatible 
with  these  bonded  phases.   Moreover,  because  of  the  wide  range  of 
polarity  of  the  mobile  phases  that  can  be  used  with  these  bonded  phases, 


an  enormous  amount  of  chemical  compounds  can  be  separated  by  using 
bonded  phases  in  RPLC. 

Despite  the  prevalence  of  RPLC  in  analytical  separations,  there  is 
a  lack  of  full  understanding  of  solute  retention  and  selectivity.   This 
lack  of  understanding  hinders  the  development  of  a  retention  index  in 
RPLC  similar  to  the  Kovats  index  (Kovats  1965)  or  Rohrschneider  and 
McReynolds  (Rohrschneider  1965;  McReynolds  1970)  constants  in  gas 
chromatography  (GC) .   It  also  creates  problems  for  practical 
chromatographers  in  developing  separation  schemes  and  for  comparison  of 
retention  results  obtained  from  different  laboratories.   Horvath  and 
coworkers  (Horvath  et  al.,  1976;  Melander  and  Horvath  1980)  explained 
their  solvophobic  retention  mechanism  for  RPLC  using  a  solute 
association  model.   This  is  an  important  step  toward  the  deconvolution 
of  solute  retention  on  the  molecular  basis  in  RPLC,  although  there  are 
many  shortcomings  in  this  theory  such  as  the  ignorance  of  the  retention 
effects  due  to  the  stationary  phase.   The  deficiencies  of  this  theory 
have  been  demonstrated  in  the  literature  by  many  researchers  (Lochmiiller 
and  Wilder,  1979;  Lochmiiller  et  al.,  1981;  Sadek  and  Carr,  1984; 
Berendsen  et  al.,  1980;  Sentell  and  Dorsey,  1989a).   Dill  and  coworkers 
(Dill  1987;  Marqusee  and  Dill,  1986a;  Dorsey  and  Dill,  1989)  have  put 
forth  a  molecular  basis  retention  mechanism  for  small  solutes  in  RPLC 
using  a  partition  model.   In  the  following  sections,  these  two 
significant  retention  mechanisms  in  RPLC  will  be  discussed  briefly. 

Solvophobic  Retention  Mechanism 
The  solvophobic  theory  was  promoted  by  Horvath  and  coworkers 
(Horvath  et  al.,  1976;  Melander  and  Horvath,  1980)  utilizing  the 


solvophobic  interaction  of  solute  association  from  Sinanoglu  (1967) .   In 
this  theory,  solute  retention  is  viewed  as  a  reversible  association  of 
the  solute  with  the  hydrocarbonaceous  ligands  of  the  stationary  phase. 
The  hydrophobic  interactions  between  the  solute  and  the  stationary  phase 
are  thought  to  come  from  the  fact  that  the  mobile  phase  in  RPLC  is 
relative  polar,  and  the  nonpolar  moiety  of  the  solute  is  repelled  and 
forced  to  associate  with  the  nonpolar  stationary  phase.   As  Melander  and 
Horvath  (1980)  stated,  the  solute  retention  in  solvophobic  theory 
actually  involves  two  steps.   The  first  step  is  the  opening  of  a  cavity 
in  the  solvent  that  has  the  same  size  and  shape  as  the  solute  molecule. 
The  second  step  is  the  placement  of  the  solute  in  the  solvent  cavity  and 
all  the  interactions  between  the  encircled  solute  and  the  solvent 
molecules  are  followed.   These  interactions  are  originated  from  the  van 
der  Waals  forces  and  the  electrostatic  interactions  amid  the  solute  and 
the  solvent  molecules. 

The  free  energy  associated  with  the  first  step  of  the  solvophobic 
theory  can  be  expressed  as 

AGe  =  KS(r)  A  y  N  (1-1) 

where  N  is  Avogadro's  number,  y   is  the  surface  tension  of  the  solvent,  A 
is  the  molecular  surface  area  of  the  solute  and  K^(r)  is  a  proportional 
factor  for  the  cavity  size.    From  equation  1-1,  the  free  energy  of  the 
solvent  cavity  is  proportional  to  the  surface  tension  of  the  solvent. 
The  free  energy  related  to  solute-solvent  interactions  is  comprised  of 
two  chemical  effects  and  an  entropic  term 


A  Gint  =  A  Gvdw  +  A  Ges  +  RT  ln(RT/PV)  (1-2) 

The  change  in  free  energy  from  the  first  chemical  effect,  A  Gvdw, 
can  be  calculated  by  using  van  der  Waals  potential  in  condensed  media 
for  the  solute  and  solvent.   The  change  in  free  energy  from  the  second 
chemical  effect,  A  Ges,  can  be  further  divided  into  dipole  interactions 
and  ionic  interactions.   The  dipole  interactions  can  be  treated 
according  to  the  Onsager  reaction  field  approach  (1936)  while  the  ionic 
interactions  can  be  estimated  using  the  Debye-Huckel  theory.   The 
entropic  term  in  equation  1-2  is  related  to  the  "free  volume"  of  the 
solute  and  this  free  volume  is  the  measure  of  the  volume  of  a  molecule 
before  it  collides  with  another  molecule.   The  solute's  molar  volume  can 
be  used  to  compute  the  free  volume  (Melander  and  Horvath  1980)  . 

As  far  as  the  stationary  phase  effects  on  the  retention  of 
solutes,  Melander  and  Horvath  (1980)  considered  them  negligible  because 
the  stationary  phase  is  nonpolar  and  the  only  attraction  between  the 
solute  and  the  stationary  phase  is  van  der  Waals  forces.   This  van  der 
Waals  force  is  insignificant  compared  to  the  van  der  Waals  interactions 
between  the  solute  and  solvent  molecules.   Although  Melander  and  Horvath 
(1980)  acknowledged  that  an  entropic  term  can  be  used  to  sum  up  the 
restricted  translational  freedom  of  the  bonded  hydrocarbonaceous  ligands 
at  the  silica  surface,  they  elected  to  pay  no  attention  to  this  term  in 
their  theory.   In  summary,  the  solvophobic  retention  mechanism  takes  the 
approach  that  solute  retention  in  RPLC  is  largely  due  to  the  hydrophobic 
interactions  of  the  solute  and  solvent  molecules.   The  stationary  phase 


is  treated  as  a  passive  entity  that  is  forced  to  receive  the  solute  from 
the  solvent,  ana  its  contribution  to  solute  retention  is  negligible. 

Partition  Retention  Mechanism 
One  of  the  catastrophic  downfalls  of  the  solvophobic  theory  is 
that  it  disregarded  the  effect  of  the  stationary  phase  on  solute 
retention.   Dill  and  coworkers  (Dill  1987;  Marqusee  and  Dill,  1986a, 
Dorsey  and  Dill,  1989)  applied  the  mean-field  statistical  mechanical 
theory,  lattice  theories  (Hill  1986)  and  random-mixing  approximation  to 
explore  the  retention  mechanism  of  RPLC  on  the  molecular  level.   Dill 
and  Dorsey  (1989)  used  a  simple  partition  model  and  took  into  account 
the  stationary  phase  effect  on  solute  retention.   They  assumed  the 
stationary  phase  in  RPLC  as  an  "interphase"  and  solute  retention  is  due 
to  partition  between  the  the  bulk  mobile  phase  and  the  interphase.   In 
their  model,  the  dominating  driving  force  of  solute  transfer  is  from  the 
differences  amid  the  chemical  affinity  of  the  mobile  and  stationary 
phase.   The  solute  capacity  factor  measured  in  RPLC  can  be  expressed  as 

k'  =  K  O  (1-3) 

where  k'  is  the  capacity  factor  of  the  solute,  K  is  the  equilibrium 
constant  of  the  partition  process  and  O  is  the  phase  ratio  of  the 
chromatographic  system  and  in  RPLC,  this  is  defined  as  the  ratio  of  the 
volumes  of  stationary  and  mobile  phases.   The  equilibrium  constant,  K, 
can  be  related  to  the  free  energy  G(S)  for  the  solute  S 

In  K  =  -(  G°sta(S)  -  G°mobile(S)  )/kT  (1-4) 


where  G°sta(S)  and  G°mo^j_le  (S)  are  the  standard-state  free  energies  of 
solute  S  in  the  stationary  phase  and  the  mobile  phase  respectively,  and 
kT  is  the  Boltzmann's  constant  multiplied  by  absolute  temperature. 

The  actual  solute  transfer  of  this  partition  process  is  comprised 
of  three  steps.   First,  a  cavity  having  the  same  size  as  the  solute 
molecule  is  opened  in  the  stationary  phase.   Second,  the  solute  molecule 
from  the  mobile  phase  is  transferred  into  the  stationary  phase  cavity. 
Third,  the  cavity  left  behind  by  the  solute  molecule  in  the  mobile  phase 
is  filled  by  other  mobile  phase  molecules.   These  three  steps  are  shown 
in  Figure  1-1,  and  for  simplification,  the  mobile  phase  is  considered  to 
be  a  single  component.   This  will  be  used  throughout  the  following 
discussion  of  this  retention  mechanism.   For  a  detailed  derivation  of 
this  retention  mechanism  using  a  binary  mobile  phase,  the  readers  should 
refer  to  Dill  (1987)  . 

The  free  energies  involved  in  these  three  steps  can  be  calculated 
using  pair  interactions  of  molecules.   The  pair  potential,  u(r*),  of 
bringing  a  spherical  molecule  X  from  an  infinite  separation  to  within  an 
average  equilibrium  separation  of  spherical  molecule  Y  can  be  defined  as 

u(r*)  -  u(°°)  =  u(r*)  =  wXY  (1-5) 

where  wXY  is  the  reversible  work.   Figure  1-2  shows  a  pair  potential, 
with  short-ranged  repulsion  and  longer-ranged  attraction.   The  partition 
process  in  RPLC  is  the  result  of  this  pair  potential.   The  exact  nature 
of  the  pair  potential  is  difficult  to  identify  since  the  partition 
coefficient  of  the  process  alone  does  not  give  enough  information,  but 
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Figure  1-1.  The  three  steps  involved  in  the  partition  process 

where  the  transfer  of  solute  molecule  S  requires  the 
opening  of  a  cavity  in  the  stationary  phase  C  and  the 
closing  of  a  cavity  in  mobile  phase  A  (Dorsey  and 
Dill,  1989)  . 


Uxv(r 


WxY  =  UxY(r*) 


Figure  1-2.  Pair  interaction  potential,  u\-v(r),  for  two  simple 

molecules.  Reversible  work  for  bringing  molecules  X 
and  Y  together  to  their  equilibrium  separation  r*  is 
wxy  (Dorsey  and  Dill,  1939) . 


it  is  generally  believed  that  it  is  related  to  hydrogen  bonding,  dipole, 
ionic  and  van  der  Waals  interactions. 

Dill  and  coworkers  (Dill  1987;  Marqusee  and  Dill,  1986a;  Dorsey 
and  Dill,  1989)  stated  that  from  the  lattice  theory  of  liquids,  the 
opening  of  the  cavity  in  the  stationary  phase  leaves  z  sides,  or  z/2 
contact  (pair  of  sides),  without  contact.   The  coordination  number,  z, 
is  the  number  of  neighbors  of  each  molecules,  or  of  each  lattice  site. 
This  cavity  opening  process  in  the  stationary  phase  contributes 
(-z/2)wq£  to  the  free  energy  of  the  overall  process.   The  same  argument 
can  be  applied  to  the  closing  of  the  cavity  in  the  mobile  phase,  and 
therefore  the  free  energy  contributed  is  (-z/2)wAA.   The  solute 
molecule,  S,  in  the  mobile  phase  is  surrounded  by  z  nearest-neighbors  of 
the  mobile  phase  molecules (solute  self-association  is  ignored  here). 
The  transfer  process  of  the  solute  from  the  mobile  phase  to  the 
stationary  phase  involves  the  formation  of  z  number  of  contacts  in  the 
stationary  phase  while  breaking  z  number  of  contacts  in  the  mobile 
phase.   The  free  energies  attributed  from  these  two  processes  are  zwgQ 
and  zwsa  respectively.   The  overall  free  energy  of  transfer  is  shown 
below. 

A  transfer  =  z  <  WSC  "  WSA  +  "AA^2  ~  wCc/2>  d"6* 

It  is  more  convenient  to  express  this  free  energy  of  transfer  in  terms 
of  the  binary  interaction  constant,  Xsc  and  XsA' 
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(^  Gtransfer>/kT  =  *SC  "  *SA  d"7' 

where  kT  is  Boltzmann's  constant  multiplied  by  the  absolute  temperature. 
The  binary  interaction  constant  is  defined  as 

Xxy   =  z/kT  <  WXY  "  (wXX  +  WYY> /2    '  (1-8) 

Equation  1-6  can  be  further  simplified  by  using  the  Hilderbrand 
solubility  parameter  concept  (Hilderbrand  and  Scott,  1950)  for  small 
molecules.   The  major  interactions  for  small  molecules  are  the  induced 
dipoles  and  the  binary  interaction  constant  can  be  approximated  as  a 
product  of  factors  involving  unitary  interaction  constants,  5X  and  8Y, 
the  solubility  parameters 

X  XY  =  constant  (  5^  -  8Y  )  (1-9) 

The  solubility  parameter  of  a  compound  can  be  measured  from  the  enthalpy 
of  vaporization  of  the  pure  compound.   Although  solubility  parameters 
have  been  used  to  model  chromatographic  retention  (Schoenmakers  et  al., 
1978;  1982;  1983;  Tijssen  et  al, .  1976;  Karger  et  al.,  1978)  and  the 
factorization  can  reduce  the  size  of  the  data  base  of  constants,  this 
approximation  is  generally  poor  compared  to  the  binary  interaction 
constant  except  in  the  case  when  dispersion  forces  are  the  principle 
interactions  (Schoenmakers  et  al.,  1978). 

The  underlying  molecular  interactions  that  drive  this  partition 
process  in  RPLC  are  very  complicated.   Since  the  capacity  factor  of  a 
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solute  in  RPLC  only  gives  information  on  the  equilibrium  constant  of  the 
partition  process,  other  experimental  data  are  needed  to  fully 
understand  these  interactions  (Dorsey  and  Dill,  1989)  .   For  example,  the 
entropic  dependence  of  these  interactions  can  be  determined  by  measuring 
the  relationship  between  temperature  and  the  capacity  factor  of  the 
solute.   The  variation  of  the  capacity  factor  with  pH  and  salt 
concentration  in  the  mobile  phase  can  be  used  to  uncover  the  effects  of 
the  electrostatic  interactions  on  the  partition  process. 

The  partition  model  is  supported  by  many  experimental  findings. 
First,  the  capacity  factors  of  solutes  in  RPLC  are  found  to  be  directly 
proportional  to  the  water/octanol  partition  coefficients  (Schantz  and 
Martire,  1987;  Opperhuizen  et  al.,  1987;  Braumann  1986;  Minick  et  al., 
1987;  Kaliszan  1986)  .   Second,  inasmuch  as  the  cavity  formation  and 
closing  in  the  stationary  and  mobile  phases  are  important,  the  capacity 
factors  of  solutes  in  RPLC  are  linearly  dependent  on  the  size  of  the 
solute  (Colin  et  al.,  1983;  Mockel  et  al.,  1987;  Horvath  et  al . ,  1976). 
Third,  the  partition  of  solutes  into  the  stationary  phase  increases  when 
the  solubility  of  the  solutes  decreases  in  the  mobile  phase  and 
therefore  the  retention  of  the  solutes  increases.   This  is  evidenced 
when  salt  content  in  the  mobile  phase  is  found  to  be  directly 
proportional  to  the  retention  of  hydrocarbons  in  RPLC  (Tanford  1980; 
Melander  and  Horvath  1980)  .   Fourth,  the  surface  tension,  yA,  of  a  pure 
mobile  phase  can  be  expressed  as 

YA  -  wAA/2a  d"10' 
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where  a  is  the  area  per  AA  contact.   Although  it  only  describes  the 
cavity  in  the  mobile  phase,  at  situations  when  the  stationary  phase 
effects  on  the  solute  retention  are  small  (i.e.,  the  other  terms  in 
equation  1-6  are  small)  the  surface  tension  of  the  mobile  phase  should 
be  proportional  to  the  retention  (Melander  and  Horvath,  1980;  Hammer  et 
al.,  1982;  Horvath  et  al.,  1976). 

Perhaps  another  possibility  of  solute  retention  in  RPLC  is  by 
adsorption.   An  adsorption  mechanism  related  to  other  types  of 
chromatography  has  been  developed  by  Martire  and  Boehm  (1983)  and 
Jaroniec  (1984) .   Dill  (1987)  has  examined  the  possibility  of  adsorption 
mechanism  and  compared  it  with  the  partition  theory  that  he  proposed. 
He  concluded  that  the  adsorption  model  is  inferior  to  the  partition 
model  for  two  reasons.   First,  in  the  adsorption  model,  the  density  of 
the  grafted  alkyl  chains  on  the  stationary  phase  should  have  no  effect 
on  the  retention  of  solutes.   On  the  other  hand,  in  the  partition  model, 
the  partitioning  of  solutes  should  decrease  if  the  density  of  the 
grafted  alkyl  chains  increases  because  it  would  be  entropically 
expensive  for  the  grafted  chain  to  uptake  the  solutes  after  a  critical 
bonding  density  (Dill  1987;  Dorsey  and  Dill,  1989)  .   This  decrease  in 
solute  partitioning  with  increase  in  alkyl  chain  density  has  been 
reported  (Claudy  et  al.,  1985;  Sentell  and  Dorsey,  1989a).   Second,  if 
partition  is  the  predominant  mechanism  of  solute  retention,  the  In  k's 
of  a  homologous  series  should  be  a  linear  function  of  the  logarithm  of 
the  appropriate  water/octanol  partition  coefficient,  with  a  slope  of  1. 
For  the  adsorption  mechanism,  such  a  plot  should  have  a  slope  of  1/z. 
This  is  because  in  adsorption,  only  one  side  of  the  solute  is  in  contact 
with  the  alkyl  chains,  and  therefore  the  transfer  free  energy  of  this 
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process  is  only  1/z  amount  of  the  partition  process  measured  by  the 
water/octanol  partition  coefficient.   In  the  partition  model,  the  solute 
is  embedded  by  the  alkyl  chains,  and  the  whole  solute  is  in  contact  with 
the  alkyl  chains.   This  simulates  the  partition  process  between  two  bulk 
liquid  phases,  and  hence  the  In  k'  of  a  solute  in  RPLC  has  a  one-to-one 
relationship  with  the  corresponding  logarithm  of  the  water/octanol 
partition  coefficient,  and  this  has  been  found  in  many  systems  (Butte  et 
al.,  1981;  Braumann  1986;  Kaliszan  1986). 

Although  Dill  and  coworkers  (Dill  1987;  Marqusee  and  Dill,  1986a; 
Dorsey  and  Dill,  1989)  used  a  few  assumptions  and  some  simple 
approximations  on  the  partition  mechanism,  their  model  correctly 
accounted  for  most  of  the  experimental  observations  that  cannot  be 
explained  sufficiently  by  other  retention  mechanisms.   This  model  does 
not  treat  the  stationary  phase  as  an  inert  part  of  the  chromatographic 
separation.   The  stationary  phase  is  considered  to  play  an  active  role 
in  the  solute  partition  process  since  the  free  energy  involved  in  the 
cavity  opening  procedure  in  the  stationary  phase  is  of  significant 
influence  on  the  total  energetic  of  the  chromatographic  separation  in 
RPLC. 

Solute-Solvent  Interaction  Constant 
The  binary  interaction  constant  %   is  an  important  parameter  for 
the  partition  mechanism  and  can  be  used  to  calculate  the  solute-solvent 
free  energy  as  seen  in  equation  1-7.   The  binary  interaction  constant  X 
is  usually  obtained  by  measuring  the  free  energy  of  transfer  or  vapor 
pressure  in  the  Henry's  law  region  (Hill  1986).   Dill  (1987)  used  the 
partition  model  and  predicted  that  in  an  RPLC  system  having  a  binary 
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mobile  phase,  the  capacity  factor  of  a  solute  can  be  estimated  by  using 
the  binary  interaction  constants  among  the  solutes  and  solvents.   He 
calculated  the  capacity  factors  of  some  solutes  by  using  the  X  obtained 
in  vapor  pressure  measurement  for  the  binary  solvents  and  X   obtained 
from  free  energy  of  transfer  of  the  solutes  and  found  good  agreement 
with  experimental  data.   The  scenario  can  be  reversed  by  using 
chromatographic  capacity  factors  of  solutes  to  evaluate  the 
solute-solvent  interaction  constant .   This  should  be  very  beneficial 
because  chromatographic  data  are  more  convenient  to  collect  and  reliable 
than  vapor  pressure  measurement  data. 

Stationary  Phase  in  RPLC 

The  retention  mechanism  in  RPLC  has  been  thoroughly  reviewed  and 
it  has  been  pointed  out  by  the  partition  mechanism  that  the  stationary 
phase  in  RPLC  contributes  a  great  deal  to  solute  retention  and  cannot  be 
ignored.   In  order  to  fully  understand  the  chromatographic  process  and 
to  forge  ahead  the  applicability  of  RPLC,  the  structure  and  properties 
of  stationary  phases  used  in  RPLC  deserve  an  in-depth  investigation. 
The  stationary  phase  in  RPLC  has  not  been  under  the  intense  study  as  the 
mobile  phase  has,  but  recently  a  few  articles  have  been  published  on  the 
structure  and  properties  of  them  (Sander  et  al.,  1983;  Gilpin  1984;  Carr 
and  Harris,  1986;  Buszewski  and  Suprynowicz,  1987;  Nawrocki  and 
Buszewski,  1988)  . 

Silica  is  by  far  the  most  common  base  material  used  in  RPLC 
stationary  phases  because  silica  is  the  best-known  inorganic,  polymeric 
material  (Unger  1979)  .   Many  reliable  methods  are  available  to  control 
the  size  and  pore  size  of  silica,  and  silica  of  the  required  size  and 
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shape  can  be  produce  in  large  quantity.   Before  the  introduction  of 
chemically  bonded  phases,  silica  was  used  as  the  support  for 
liquid-liquid  chromatography  where  a  thin  film  of  liquid  is  coated  onto 
the  surface  of  the  silica  and  then  is  used  as  the  stationary  phase 
(Unger  1979)  .   These  liquid  stationary  phases  suffer  from  instability 
and  irreproducibility .   The  evolvement  of  chemically  bonded  phases  by 
chemically  treating  the  silica  surface  with  silanes  began  in  the  middle 
of  1960s.   They  were  debuted  in  gas  chromatography  as  nonpolar 
stationary  phase  obtained  by  reacting  silica  gel  with 
hexadecyltrichlorosilane  (Abel  et  al.,  1966).   The  employment  of 
chemically  bonded  phases  in  liquid  chromatography  was  first  suggested  by 
Stewart  and  Perry  (1968) .   Chemically  bonded  phases  have  excelled  to 
become  the  most  frequently  used  stationary  phases  in  RPLC  today  due  to 
the  advancement  of  synthetic  technology  for  preparing  these  phases. 

The  underlying  properties  of  these  chemically  bonded  phases  depend 
on  the  physical  and  chemical  properties  of  the  base  silica,  the  bonding 
reactions  and  reagents  used.   The  base  silica  constitutes  the  major 
portion  of  the  stationary  phase,  and  the  first  property  of  the  base 
silica  that  influences  the  chromatographic  process  is  the  particle  size 
of  the  silica.   The  size  of  most  silica  particles  employed  in  RPLC 
stationary  phases  are  in  the  range  of  3-10  (im.   The  actual  effect  of  the 
silica  particle  size  on  RPLC  is  inconclusive  since  inconsistent  results 
have  been  reported  in  the  literature.   Dewaele  and  Verzele  (1983)  have 
done  a  study  on  the  influence  of  silica  particle  size  distribution  on 
the  reversed  phase  packings,  and  they  found  that  the  particle  size 
affects  chromatographic  performance  of  the  packings  minimally.   Gazda  et 
al.  (1980)  also  had  the  same  conclusions,  but  Bristow  (1978)  has  shown 
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otherwise.   Pore  size  distribution  is  another  important  physical 
property  of  silica  particles.   Verzele  et  al.  (1985)  have  noticed  that 
the  mean  pore  diameter  of  the  silica  used  as  chromatographic  packings 
should  be  in  the  range  of  6-10  nm,  and  the  narrower  the  distribution 
curve  the  better  it  is.   Micropores  of  diameter  2  nm  or  smaller  are 
undesirable  because  they  can  alter  the  final  coverage  of  the  bonded 
phase  (Berendsen  et  al.,  1980;  Engelhart  et  al.,  1982). 

Apart  from  the  physical  properties,  the  surface  chemistry  of  the 
silica  particles  endure  the  final  characteristic  of  the  bonded  phase. 
Silica  surface  is  very  complex,  and  it  has  been  found  to  consist  of 
various  kinds  of  silanols  and  siloxane  bonds.   The  siloxane  bonds  are 
usually  considered  as  rather  chemically  inert  and  therefore  do  not  take 
part  in  the  formation  of  the  bonded  phase.   Single,  geminal  and  vicinal 
forms  of  silanols  are  found  on  the  surface  of  silica  (Unger  1979;  Majors 
1980) ,  and  chromatographers  have  mixed  feelings  toward  these  surface 
silanols.   First  of  all,  these  silanols  react  with  silanes  to  form  the 
bonded  phase,  but  on  the  other  hand,  these  silanols  are  one  of  the  major 
contributions  to  the  variability  of  the  bonded  phases.  Moreover,  these 
silanols  account  for  the  pH  of  the  silica.   Different  values  of  pKa  of 
silica  have  been  reported  in  the  literature  (Walling  1950;  Karger  et 
al.,  1980;  Majors  1980).   Engelhardt  and  Muller  (1981)  first  observed 
the  divergence  of  surface  pH  of  silica  packings  for  chromatographic 
purpose.   This  reflects  the  different  methods  of  manufacturing  silicas 
and  connects  to  the  different  chromatographic  properties  of  silicas 
(Hansen  et  al.,  1986)  . 

Trace  metal  impurities  found  in  silica  is  another  factor  which 
influences  the  property  of  the  bonded  phase.   The  effect  of  trace  metal 
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impurities  have  been  ignored  by  many  researchers,  but  Marshall  et  al. 
(1986;  1987),  Nawrocki  (1986;  1987)  and  Sadek  et  al.  (1987)  have  shown 
that  these  impurities  can  indeed  affect  the  chromatographic  process. 
Verzele  et  al.  (1979)  have  pointed  out  that  numerous  transition  metals 
can  be  found  in  the  silicas  that  are  used  for  chromatographic  packings. 
Sadek  et  al.  (1987)  described  that  silica  can  contain  metals  in  three 
forms:  surface  species,  internal  and  secluded.   The  surface  metals  can 
form  surface  metal  hydroxides  which  can  alter  the  surface  pH  of  the 
silicas,  while  both  the  surface  and  the  internal  metals  can  participate 
in  coordination  with  the  solute  or  solvents  using  their  orbitals.   The 
secluded  metals  are  too  far  away  from  the  surface  silanols  and  therefore 
do  not  exert  any  effect  on  the  chromatographic  process .   According  to 
Verzele  and  coworkers  (Verzele  et  al.,  1979;  Verzele  and  Dewaele,  1981; 
Verzele  1983)  and  Sadek  et  al.  (1987),  by  acid  treatment  of  the  silica, 
most  of  the  metal  impurities  can  be  leached. 

Chemically  bonded  phases  are  silica  particles  having  their  surface 
derivatized  with  alkyl  chains  or  other  functionalities  such  as  cyano  or 
amino  groups  (Dorsey  and  Dill,  1989) .   The  surface  of  the  silica 
particles  can  be  derivatized  by  reacting  with  organochlorosilanes  or 
organoalkyoxysilanes,  and  the  most  common  synthetic  scheme  is  called 
"monomeric  reaction".   In  this  type  of  reaction,  a  single  siloxane 
linkage  is  formed  between  one  silanol  on  the  surface  of  the  silica  and  a 
monochlorosilane  or  monoalkoxysilane .   A  simplified  diagram  showing  this 
type  of  monomeric  reaction  is  depicted  in  Figure  1-3.   The  popularity  of 
monomeric  reaction  is  due  to  the  fact  that  it  gives  a  uniform  and 
well-defined  single  layer  coverage  of  the  surface  of  the  silica,  and 
also  the  reproducibility  is  above  other  types  of  reactions.   The  maximum 
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number  of  silanols  on  the  surface  of  silica  are  predicted  to  be  about  8 
Jlmol/m^  and  only  these  silanols  can  react  with  the  silanes  (Unger  1979). 
The  most  common  Cg  and  C^g  containing  silanes  usually  produce  a  surface 
coverage  of  approximately  2.5-3.0  fimol/m  .   This  means  that  there  is  an 
immense  amount  of  residual  silanols  which  does  not  react  with  the 
silanes.   These  residual  silanols  are  accessible  to  solutes  during  the 
chromatographic  process  and  can  cause  adsorption  and  mixed  retention 
mechanism.   A  second  bonding  procedure  is  usually  employed  by  most 
manufacturers  to  try  to  minimize  the  residual  silanols.   The  second 
bonding  procedure  often  makes  use  of  trimethylsilanes  to  react  with  the 
remaining  silanols  and  is  termed  "end-capping"  (Lochmuller  and  Marshall, 
1982;  Sadek  and  Carr,  1983;  Evans  et  al.,  1980;  Dewaele  et  al.,  1982; 
Marshall  et  al.,  1987).   Other  methods  include  using  silica  pretreatment 
(Gobet  and  Kovats,  1984;  Marshall  et  al.,  1984),  and  the  use  of  more 
reactive  silanes  have  been  reported  (Buszewski  et  al.,  1987;  Lork  et 
al.,  1986) . 

Much  research  has  been  done  on  improving  the  synthetic  means  to 
increase  the  surface  coverage  or  bonding  density  of  the  monomeric  bonded 
phase  (Buszewski  et  al.,  1987;  Golding  et  al.,  1987;  Khong  and  Simpson, 
1987;  Sentell  et  al.,  1988) .   This  striving  for  higher  bonding  density 
originated  as  a  way  to  suppress  the  residual  silanols,  but  it  was  soon 
noticed  that  there  are  other  advantages  to  high  bonding  density.   First, 
by  increasing  the  bonding  density,  less  silanols  are  exposed  to  the 
solutes  and  thus  eliminates  most  of  the  adsorption  mechanism  due  to  the 
silanols.   Second,  the  stability  of  the  stationary  phase  can  be  raised 
by  higher  bonding  density  because  the  dense  alkyl  chains  protect  the 
silica  from  being  hydrolyzed. 
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When  a  di-  or  trifunctional  silane  is  used  instead  of  a 
monofunctional  silane  in  the  bonding  process,  a  polymeric  coverage 
frequently  results.   The  polymeric  network  is  a  consequence  of  the 
reaction  among  one  or  more  of  the  leaving  groups  of  the  silane  and  the 
silanols  on  the  silica  surface.   Sander  and  Wise  (1984)  have  noted  that 
under  a  controlled  manner,  polymeric  phases  can  be  very  reproducible  and 
have  some  chromatographic  advantages.   Wise  and  Sander  (1985)  showed 
that  these  polymeric  phases  have  shape  selectivity  on  polycyclic 
aromatic  hydrocarbons  and  are  similar  to  liquid  crystalline  phases  used 
in  GC .   They  proposed  that  polymeric  phases  are  more  ordered  than 
monomeric  phases. 

The  microscopic  structure  of  the  bonded  phase  is  very  important  to 
the  retention  and  selectivity  behavior  of  the  bonded  phase.   Many 
research  groups  have  tried  many  different  experimental  approaches  to 
gain  better  understanding  of  the  nature  and  molecular  organization  of 
the  alkyl  chains  that  have  been  grafted  onto  the  silica.   In  the  early 
stage  of  the  bonded  phase,  a  simplified  model  was  suggested  and  is  shown 
in  Figure  1-4  (Melander  and  Horvath,  1980) .   This  model  is  not  very 
realistic  because  it  assumes  that  the  alkyl  chains  are  stiff  rods 
without  any  internal  degrees  of  freedom.   This  is  contradictory  to  the 
corresponding  disorder  structure  of  alkane  at  the  temperature  of 
interest  for  chromatography.   Also,  the  "fur"  and  "stack"  models  in 
Figures  l-4b  and  l-4c  imply  that  the  alkyl  chains  are  exposed  to  the 
mobile  phase,  but  the  mobile  phase  used  in  RPLC  is  of  very  high  aqueous 
content  and  the  hydrophobic  effect  should  prohibit  this  total  exposition 
of  the  nonpolar  alkyl  chains  to  the  polar  mobile  phase  (Dill  1987)  . 
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Early  models  of  molecular  structure  and  organization 
of  the  bonded  phase  in  RPLC:  a)  "picket  fence"; 
b)  "fur";  c)  "stacks"  (Dill  1987)  . 
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Figure    1-5.    Interphase   model    of    the   bonded   phase    in    RPLC   proposed 
bv   Dill    (1987) . 
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Dill  and  coworkers  (Dill  1987;  Marqusee  and  Dill,  1986a;  Dorsey 
and  Dill,  1989)  have  proposed  an  interphase  model  for  the  structure  of 
the  bonded  phase.   Figure  1-5  is  a  simple  representative  of  this  model. 
The  alkyl  chains  that  are  bonded  onto  the  surface  should  have  the  same 
molecular  structure  and  organization  of  other  similar  interfacial  phases 
involving  chain  molecules  such  as  micelles,  surfactant  monolayers, 
bilayers  and  microemulsions  (Dill  and  Flory,  1980;  1981;  Dill  et  al., 
1984) .   A  typical  interphase  has  one  side  of  the  chain  molecules 
anchored  to  an  interface,  and  the  disorder  of  the  chains  increases  with 
the  distance  from  the  interface  (Dill  et  al.,  1988;  Brown  1984;  Cabane 
1977) .   The  bonded  phase  in  RPLC  has  these  characteristics  because  one 
side  of  it  is  the  siloxane  bond  between  the  alkyl  chains  and  the  surface 
silanol.   Gilpin  and  coworkers  (Gilpin  and  Gangoda,  1983;  1984;  Gilpin 
1984)  have  shown  by  using  NMR  measurements  of  the  T]_  of  1-^C  labelled 
alkyl  chains  that  segmental  motion  of  the  alkyl  chains  is  a  function  of 
distance  from  point  of  surface  attachment.   Sander  et  al.  (1983)  used 
fourier  transform  infrared  spectroscopy  to  study  a  series  of 
dimethyl-n-alkyl  bonded  phases  ranging  from  C]_  to  C22 •   Their  results 

showed  that  these  systems  have  disordered  chains  with  kinks  and  bends. 

The  molecular  organization  of  the  interphase  is  dominated  by  three 
factors  (Dill  and  Flory,  1980;  1981;  Dill  et  al.,  1984;  Marqusee  and 
Dill,  1986b) .   First,  the  molecular  organization  of  the  interphase  is 
restricted  by  the  geometry  of  the  interface  and  the  density  and  length 
of  the  chain  molecules.   Second,  under  a  poor  wetting  agent,  the 
interphase  will  try  to  reject  most  of  the  hostile  solvent  molecules.   An 
example  is  when  alkyl  stationary  phase  of  RPLC  is  under  aqueous  mobile 


24 

phase.   Third,  in  accordance  with  the  attempt  to  gain  maximum  entropy, 
the  interphase  always  adopts  as  much  disorder  as  possible  within  all  the 
constraints . 

The  interphase  model  of  the  stationary  phase  is  very  appropriate 
with  the  partition  mechanism  of  solute  retention.   In  theory,  the 
partition  of  a  solute  molecule  into  the  interphase  having  a  fixed 
surface  density  will  cause  ordering  of  the  alkyl  chains  (Margusee  and 
Dill,  1986a;  Dill  1987;  Dorsey  and  Dill,  1989),  and  therefore  solute 
retention  is  entropically  unfavorable.   The  partition  of  the  solute 
should  be  directly  proportional  to  the  surface  coverage  of  the  alkyl 
chains  until  it  reaches  a  maximum  where  interactions  between  individual 
chains  become  significant.   Then  the  partition  coefficient  will  decrease 
to  zero  when  the  maximum  chain  density  is  reached,  which  is  roughly  8 
umol/m  .   The  work  by  Sentell  and  Dorsey  (1989a)  using  bonded  phases 
having  different  bonding  densities  has  confirmed  this  prediction.   Since 
the  end  of  the  alkyl  chains  on  the  bonded  phase  are  found  to  be  in  rapid 
motion  compared  to  their  attached  ends,  the  solute  will  prefer  to 
partition  to  the  free  ends  of  the  chains.   White  et  al.  (1981)  observed 
this  type  of  preferential  distribution  when  they  performed  neutron 
scattering  experiments  on  similar  bilayer  membranes. 

The  importance  of  the  interphase  model  is  that  it  predicts  the 
essential  of  the  stationary  phase  effects  on  solute  retention  and 
selectivity  in  RPLC.   Dependence  of  solute  retention  on  the  surface 
density  of  alkyl  chains  can  be  explained  by  the  model.   Selectivity  due 
to  different  molecular  shapes  have  been  widely  noticed  (Wise  and  Sander, 
1985;  Tanaka  et  al.,  1980;  Tanaka  et  al.,  1982;  Lochmuller  et  al., 
1985) .   Using  the  interphase  model,  this  shape  selectivity  can  be 
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considered  as  arising  from  the  fact  that  more  free  energy  is  required  to 
align  solutes  having  shapes  parallel  to  the  alkyl  chains  compared  to 
molecules  that  have  shapes  normal  to  the  alkyl  chains.   Sentell  and 
Dorsey  (1989b)  measured  selectivity  of  six  four-ring  polycyclic  aromatic 
hydrocarbons  (PAH)  on  stationary  phases  having  bonding  densities  from 
1.74-4.07  ^mol/m2.   They  found  that  the  selectivity  of  the  PAHs  increase 
swith  the  surface  bonding  density.   The  increase  in  bonding  density 
causes  ordering  of  the  alkyl  chains  and  hence  elevates  the  shape 
selectivity  of  these  PAHs. 

In  this  study,  the  partition  mechanism  proposed  by  Dill  and 
coworkers  (Dill  1987;  Marquee  and  Dill,  1986a;  Dorsey  and  Dill,  1989)  is 
tested  against  an  extended  data  base.   Results  will  be  presented  in 
Chapter  II  to  show  that  the  partition  mechanism  is  correct  and  the 
solute-solvent  free  energy  of  the  binary  mobile  phase  and  the  solutes 
can  be  obtained  from  chromatographic  data  as  predicted  from  the  theory. 
The  ample  number  of  stationary  phases  for  RPLC  available  today  all  have 
distinct  retention  and  selectivity  properties.   This  is  likely  exerted 
by  the  different  stationary  phase  structure  and  properties  on  the  solute 
retention  mechanism.   Many  different  stationary  phases  were 
investigated,  and  the  relative  retention  strength  of  them  were  found  to 
be  related  to  the  phase  ratio  <t>,    and  the  functionality  of  the  silane 
used  in  forming  the  stationary  phase.   Results  of  this  study  will  be 
introduced  in  Chapter  III  along  with  a  simple  and  useful  method  that  was 
derived  to  classify  these  stationary  phases  according  to  their 
retentivity  towards  all  solutes. 


CHAPTER  II 
SOLUTE-SOLVENT  INTERACTION  FREE  ENERGIES 
IN  REVERSED  PHASE  LIQUID  CHROMATOGRAPHY 

Introduction 

Binary  mobile  phases  of  water  and  an  organic  modifier  such  as 

acetonitrile,  methanol  or  THF  are  most  commonly  used  in  RPLC.   From  the 

partition  model  proposed  by  Dill  and  coworkers  (Dill  1987;  Marqusee  and 

Dill,  1986a,  Dorsey  and  Dill,  1989) ,  the  equilibrium  constant,  K,  of  a 

binary  mobile  phase  with  components  A  and  B  (a  case  not  discussed  in 

detail  in  the  previous  chapter)  can  be  defined  as  a  simple  quadratic 

function  of  the  mobile  phase  composition  and  the  binary  interaction 

constants  between  the  mobile  phase  components  and  the  solute 

ln  KC/AB  "  <XSA  "  *SC>  +  <t>B  <XSB  "  *SA  "  *AB>  +  <t>B2  *AB    J2"1) 

where  0  <,   <))B  <  1  is  the  fraction  of  the  mobile  phase  sites  occupied  by  B 
molecules,  C  is  the  interphase  stationary  phase,  XsA'  ZSB'  anc*  Xsc  are 
the  binary  interaction  constants  between  the  solute  and  the  components  A 
(water) ,  B  (organic  modifier) ,  and  C  (stationary  phase)  of  the 
chromatographic  system  respectively,  and  Xab  is  the  binary  interaction 
constant  between  the  mobile  phase  components  A  and  B. 

The  quadratic  relationship  between  the  equilibrium  constant  and 
the  volume  fraction  of  organic  modifier  has  been  previously  reported  by 
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Schoenmakers  and  his  colleagues  (Schoenmakers  et  al.,  1978;  1982;  1983) 
using  the  solubility  parameter  theory  (Hilderbrand  and  Scott, 1950). 
Dill  and  coworkers  (Dill  1987;  Dorsey  and  Dill,  1989)  have  suggested 
that  instead  of  the  quadratic  equation,  a  more  useful  linear  expression 
can  be  employed  to  plot  the  dependence  of  retention  on  mobile  phase 
composition.   Since  k'/kw  =  K/Kw  where  the  subscript  w  refers  to  the 
mobile  phase  composition  when  4>b=0,  equation  2-1  can  be  rearranged  to 

1/({)B  In  (k'/k„)  =  (XSB  "  ZsA  -  *AB>  +  ^B  *AB  <2~2> 

As  stated  in  Chapter  I,  the  binary  interaction  constant  %   is  often 
determined  by  measuring  the  free  energy  of  transfer  or  vapor  pressure  in 
the  Henry's  law  region  (Hill  1986).   A  careful  examination  of  equation 
2-2  reveals  that  a  plot  of  1/<(>B  In  (k'/kw)  versus  (J)B  should  be  linear  if 
all  the  assumptions  made  in  the  partition  model  hold.   Also,  the  slope 
of  such  a  plot  should  give  the  binary  interaction  constant  XaB'  anci  the 

y-intercept  at  §B=1    should  be  equal  to  (%sb  ~  XSA*  '   This  predicts  that 
the  binary  interaction  constants  of  mobile  phases  used  in  RPLC  can  be 
obtained  from  chromatographic  data,  and  by  applying  equation  1-7,  the 
solute-solvent  free  energies  between  the  solute  and  the  two  components 
of  the  mobile  phase  in  RPLC  can  be  calculated. 

In  the  original  derivation  of  the  partition  model  (Dill  1987) ,  the 
solute  was  treated  as  having  a  size  comparable  to  the  mobile  phase 
components.   In  actuality,  solutes  in  RPLC  are  frequently  larger  than 
the  mobile  phase  components.   Since  the  cavity  term  is  very  important  in 
the  partition  model,  solutes  of  different  sizes  should  create  cavities 
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of  different;  magnitudes.   If  the  solute  molecules  are  larger  than 
molecules  of  the  mobile  phase  components,  equation  2-2  can  be  redefined 


l/63  In  (k'/kw)  =  n  [ (XsB  -  ZSA  "  *AB>  +  ^B  Zab^  <2"3> 

where  n  is  the  ratio  of  the  size  of  solute  to  the  size  of  the  solvent 
molecules.   This  expression  results  because  the  number  of  contacts  of 
solute  with  the  mobile  phase  components  and  the  size  of  the  cavities  are 
directly  proportional  to  the  size  of  the  solute  molecule. 

Theoretical  derivations  are  not  valid  until  they  have  been  proven 
experimentally.   With  this  in  mind,  we  tested  the  prediction  of  the 
composition  dependence  expressed  in  equations  2-2  and  2-3  with  an 
extensive  data  base  consisting  of  more  than  300  sets  of  experimental 
retention  data  of  various  solutes  on  various  RPLC  columns.   We  hope  to 
uncover  the  underlying  physical  nature  of  the  slope  and  intercept  of 
this  type  of  composition  plot.   In  agreement  with  the  above  predictions, 
linear  dependence  of  In  k'  on  surface  area  or  van  der  Waals  volume  of 
solute  molecules,  and  hence  the  size  of  the  cavity,  has  been  widely 
observed  (Jinno  and  Kawasaki,  1983a;  1983b;  Funasaki  et  al.,  1986;  Arai 
et  al.,  1987;  Kaliszan  1986;  Mockel  et  al.,  1987). 

Experimental  Procedure 
We  tested  the  above  predictions  with  a  large  data  base  that  had 
been  previously  used  for  correlations  of  solvatochromatic  solvent 
polarity  measurements  with  reversed  phase  retention  (Johnson  et  al., 
1986) .   The  data  base  is  used  here  in  its  entirety,  with  the  omission  of 
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six  data  sets  which  had  retention  values  at  only  three  mobile  phase 
compositions.   Additional  data  sets  generated  in  our  laboratory  (Michels 
1989)  with  ethanol-water  and  propanol-water  mobile  phases  are  also 
included. 

Linear  regression  calculations  were  performed  using  the  program 
CURVE  FITTER  (Interactive  Microware,  State  College,  PA,  USA)  on  an  Apple 
II  Plus  48K  microcomputer.   The  program  was  used  to  calculate  the 
coefficients  of  correlation,  r,  and  determination,  r  ,  of  the  data. 
This  program  was  also  used  to  extrapolate  k'  to  the  composition  of  pure 
water,  denoted  kw.   The  In  kw  values  are  not  usually  reported  in  the 
literature  as  they  are  not  easily  experimentally  determined.   We 
extrapolated  plots  of  In  k1  versus  ET(30)  polarity  to  the  polarity  value 
of  100%  water  to  obtain  the  In  kw  values  used  here.   This  method  has 
been  shown  to  give  more  reliable  values  compared  to  other  types  of 
extrapolation  (Snyder  et  al.,  1979;  Dolan  et  al.,  1979;  Schoemakers  et 
al.,  1979;  Antle  et  al.,  1985;  Baty  and  Sharp,  1988)  because  the 
calculated  In  kw  for  a  given  solute  is  consistent  among  methanol-water, 
ethanol-water,  and  acetonitrile-water  mobile  phases  when  using  this 
method  (Michels  1989) .   The  F-values,  which  are  used  to  compare  the 
significance  level  of  two  variances,  were  obtained  from  the  StatWorks 
(Heydon  and  Son  Inc.,  Philadelphia,  PA,  USA)  program  on  a  Macintosh  SE 
(Apple  Computer,  Inc.,  Cupertino,  CA,  USA)  microcomputer.   The  van  der 
Waals  volume  of  the  solutes  were  calculated  using  methods  reported  by 
Bondi  (1964),  and  the  hydrocarbonaceous  surface  area  was  obtained  from 
Woodburn  (1985) . 
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Results  and  Discussion 
The  extensive  data  base  presented  here  permits  certain  tests  of 
the  cavity  contributions  to  the  partition  retention  model,  and  the 
ability  of  using  chromatographic  data  to  calculate  binary  interaction 
constants.   A  typical  example  of  this  type  of  composition  plot  is  shown 
in  Figure  2-1.   In  the  case  shown,  the  solute  is  4-nitrophenol  and  the 
mobile  phase  is  a  series  of  different  acetonitrile-water  mixtures.   From 
the  figure,  it  is  obvious  that  this  plot  provides  a  linear 
representation  of  the  experimental  data.   Table  2-1  presents  an 
extensive  compilation  of  the  results  of  plotting  data  in  this  form  for  a 
wide  range  of  solutes,  stationary  phases,  and  various  mobile  phase 
mixtures  commonly  used  in  RPLC.   From  the  correlation  coefficients 
shown,  it  is  clear  that  the  degree  of  linearity  of  this  type  of  plot  is 
extremely  good,  with  only  a  very  small  number  of  exceptions  (Figure 
2-2) .   More  than  80%  of  the  data  sets  are  found  to  have  r^  of  0.9  and 
higher  and  almost  50%  with  r^  of  0.99  and  higher.   An  F-value,  which  is 
used  to  compare  the  significance  level  of  fitting  the  data  sets  to  a 
quadratic  model  (Anderson  1987),  is  calculated  and  expressed  as  a  in 
Table  2-1.   The  results  of  the  F-test  reveal  that  the  quadratic  model 
does  not  fit  the  data  sets  significantly  better  (at  the  90%  confidence 
level),  except  in  148  out  of  346  data  sets.   We  conclude  that  the  linear 
fitting  of  the  data  is  adequate.   The  principal  exception  in  the  data 
set  is  the  Sepralyte  C-2  column  with  methanol  as  modifier.   The 
coefficient  of  determination,  r  ,  for  this  particular  data  set  varies 
from  0.4460  to  0.9162,  with  an  average  value  of  0.7287.   These  poor 
correlations  are  probably  observed  because  for  the  short-chain 
stationary  phases,  adsorption  is  expected  to  be  the  dominant  retention 
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l/63  ln(k'/kw) 


1.0 


Figure  2-1.  Plot  of  1/<})B  ln(k'/kw)  versus  <J)B  for  the  solute 
4-Nitrophenol  using  Hypersil  ODS  column  and 
acetonitrile  as  modifier. 
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Figure  2-2.  Histogram  of  coefficient  of  determination,  r2,  for  the 
plots  of  1/<}>b  ln(k'/k;w)  versus  <J>b  for  all  the  data 
sets  . 
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mechanism  rather  than  partitioning,  which  is  the  basis  of  the  assumption 
of  equations  2-2  and  2-3.   Another  possibility  is  that  the  residual 
silanol  groups  are  more  accessible  to  solutes  in  short-chain  stationary 
phases  than  for  longer  chains,  and  they  may  exert  a  strong  effect  on 
retention  that  is  not  treated  in  the  partition  model.   Due  to  the 
significant  curvature  and  the  poor  correlation  found  in  this  data  set, 
we  feel  that  the  y-intercepts  at  <t>B=l  are  unreliable,  and  therefore  are 

not  reported  in  Table  2-1.   Where  linearity  is  found  to  be  good,  it 
implies  that  the  partition  model,  and  the  few  assumptions  applied  are 
valid.   This  linearity  is  expected  to  be  limited  to  the  intermediate 
compositions  explored  here,  because  the  random-mixing  approximation 
assumed  in  the  original  derivation  is  expected  to  fail  in  the  limits  of 
extreme  mobile  phase  composition  of  a  few  percent  of  either  mobile  phase 
component  (Schoenmakers  et  al.,  1983)  .   In  those  cases,  solutes  or  the 
minor  component  of  the  mobile  phase  may  associate  to  form  non-random 
mixtures . 

According  to  equation  2-3,  the  slope  of  the  regression  of  this 
composition  plot  equals  the  binary  interaction  constant  characterizing 
the  pair  interaction  between  molecules  of  the  A  and  B  components  in  the 
mobile  phase,  multiplied  by  the  size  of  the  cavity  occupied  by  the 
solute.   Tests  of  this  prediction  of  the  cavity-size  dependence  are 
shown  in  Figures  2-3  and  2-4.   The  slopes  of  the  composition  plots  are 
found  to  increase  with  increasing  cavity  size  created  by  the  solute 
molecules.   That  is,  slope  of  1/<1>b  In  (k'/kw)  versus  (J>B  =  nX^B  ancl 

3slope  /  dn   =  XAB  (2-5) 
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equation  2-3 
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Figure  2-3.  Plot  of  slopes  from  equation  2-3  versus  van  der  Waals 
volume,  Vw,  of  the  solutes  for  the  Sepralyte  C-18 
column  with  acetonitrile  as  modifier. 
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Table  2-2.      Regression  results  of  slopes  from  equation  2-3  versus  the 
van  der  Waals  volume,  Vw,  of  the  solutes  for  all  the 
columns . 


Column  Modifier         s^lope  (xlQ2_L  r2 

Sepralyte  C-2  Acetonitrile  16.51  0.9561 

Sepralyte  C-4  Acetonitrile  16.61  0.9263 

Sepralyte  C-8  Acetonitrile  7.817  0.7897 

Sepralyte  C-18  Acetonitrile  9.582  0.9164 

Sepralyte  C-4  Methanol  7.286  0.9189 

Sepralyte  C-8  Methanol  6.228  0.9305 

Hypersil  ODS  Acetonitrile  12.78  0.5499 

Unisil  Q  C-18  Acetonitrile  8.640  0.4106 

Ultrasphere  ODS  Methanol  10.33  0.9977 

Ultrasphere  ODS  Ethanol  5.948  0.2490 

Ultrasphere  ODS  1-Propanol  4.662  0.0495 

Ultrasphere  ODS  Acetonitrile  18.57  0.8091 

Silasorb  C-8  Methanol  3.051  0.2339 

Silasorb  C-8  Acetonitrile  2.595  0.2083 
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Slopes  from 
equation  2-3 
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Figure  2-4.  Plot  of  slopes  from  equation  2-3  versus 

hydrocarbonaceous  surface  area,  HSA,  of  the  solutes 
for  the  Sepralyte  C-18  column  with  acetonitrile  as 
modifier . 
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Table  2-3.      Regression  results  of  slopes  from  equation  2-3  vs  the 

hydrocarbonaceous  surface  area,  HSA,  of  the  solutes  for 

all  the  columns . 


Column  Modifier         Slope (xlQ2!  £2 

Sepralyte  C-2  Acetonitrile  8.511  0.9642 

Sepralyte  C-4  Acetonitrile  8.785  0.9465 

Sepralyte  C-8  Acetonitrile  4.115  0.8226 

Sepralyte  C-18  Acetonitrile  5.099  0.9280 

Sepralyte  C-4  Methanol  3.817  0.9068 

Sepralyte  C-8  Methanol  3.358  0.9520 
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where  slope  is  defined  as  the  slope  of  equation  2-3.   Values  of  XaB 
generated  from  the  derivative  of  equation  2-5  are  presented  in  Tables 
2-2  and  2-3.   Due  to  the  limited  range  of  solute  sizes  available,  this 
linearity  is  observed  when  either  the  cavity  volume  (Figure  2-3;  slopes 
in  Table  2-2),  as  measured  by  the  van  der  Waals  volume  of  the  solute 
molecule  or  cavity  area,  as  measured  by  the  hydrocarbonaceous  surface 
area  of  the  solute  molecule  (Figure  2-4;  slopes  in  Table  2-3), 
represented  the  size  of  the  cavity  created  by  the  solute  molecule.   The 
binary  interaction  constants  have  different  units  depending  on  whether 
cavity  size  is  taken  to  be  that  of  the  solute  volume  or  area.   The 
linear  relationship  holds  well  between  the  slopes  from  equation  2-3  and 
the  cavity  size;  with  most  of  these  plots  having  r2  of  0.9  and  better. 
This  linearity  implies  that  the  slopes  from  equation  2-3  are 
proportional  to  the  size  of  the  cavity  occupied  by  the  solutes.   The 
slopes  from  Figures  2-2  and  2-3  represent  the  binary  interaction 
constant  of  the  mobile  phase  component  A  with  B  per  unit  volume  or  area 
respectively.   For  the  Sepralyte  C-18  column,  shown  in  Figure  2-4  and 
Table  2-3,  the  interaction  free  energy  of  the  mobile  phase  components  is 
found  to  equal 

(5.099  x  10"2)kT  =  (5.099  x  10"2) (592  cal/mol)  =  30.2  cal/mol  A2   <2_6> 

for  these  data  taken  at  25°C. 

Several  columns  give  poor  correlations  for  the  slope  of  equation 
2-3  versus  van  der  Waals  volume  of  the  solutes.   For  the  case  of  the 
Hypersil  column  with  acetonitrile  as  modifier,  most  of  the  solutes  were 
nitrogenous  compounds  such  as  aniline,  pyridine  and  amino-substituted 
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PAHs .    It  has  been  shown  that  many  compounds  with  a  nitrogenous  moeity 
often  exhibit  an  unsatisfactory  degree  of  peak  tailing  in  RPLC,  and  this 
effect  is  generally  considered  to  come  from  the  presence  of  residual 
silanol  groups  on  the  surface  of  the  bonded  phase  (Wahlund  and 
Sokolowski,  1978;  Bayer  and  Paulus,  1987;  Smith  et  al.,  1986;  Hansen  et 
al.,  1988).   The  retention  mechanism  of  these  compounds  is  most  likely  a 
mixed  partition/adsorption  process  because  of  the  interactions  with  the 
residual  silanol  groups,  and  would  not  be  expected  to  fit  the  theory 
tested  here.   The  Unisil  Q  C-18  data  include  many  isomers  of  substituted 
benzene.   The  van  der  Waals  volumes  calculated  from  Bondi's  (1964) 
method  for  geometric  isomers  are  identical,  but  they  give  different 
slopes  for  equation  2-3.   It  has  been  shown  in  the  literature  that 
isomeric  alkylbenzenes  give  relatively  poor  correlation  with  van  der 
Waals  volume  (Smith  1981,  Jinno  and  Kawasaki,  1983a) .   The  poor 
correlation  in  this  instance  most  likely  arises  from  the 

misrepresentation  of  the  calculated  van  der  Waals  volume,  rather  than  a 
breakdown  of  the  retention  theory.   The  poor  correlations  for  the 
ethanol  and  propanol  data  on  the  Ultrasphere  ODS  column  may  be  caused  by 
the  changing  stationary  phase  environment  as  these  more  hydrophobic 
modifiers  are  used  in  the  mobile  phase.   There  is  not  yet  any  theory 
published  on  the  uptake  of  solvent  by  the  stationary  phase  chains. 
Michels  (1989)  has  recently  reported  that  propanol  appears  to  saturate 
the  stationary  phase  at  a  very  low  mobile  phase  percentage  of  propanol. 
The  poor  correlations  with  the  Silasorb  column  have  yet  to  be  accounted 
for. 

The  y-intercepts  at  <J>b=1  of  equation  2-3  are  predicted  to  equal 
the  free  energy  of  transfer  of  the  solute  between  pure  A  and  pure  B  of 
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the  mobile  phase,  multiplied  by  the  cavity  size.   That  is,  the 
y-intercepts  of  1/<{»B  In  (k'/kw)  versus  (^  at  (0B=D  =  n  (Xsb  -  Xsa>  and 

3y-intercepts  (4>B=1)  /  3n  =  (XSB  -  Xsa>  (2-7) 

The  dependence  of  these  intercepts  on  the  cavity  volume  and  area 
can  be  seen  in  Figures  2-5  and  2-6  respectively.   The  correlation 
coefficients  of  these  intercepts  and  the  cavity  size  are  generally 
greater  than  0.9,  and  the  regression  results  are  listed  in  Tables  2-4 
and  2-5.   The  intercepts  are  found  to  increase  linearly  with  cavity  size 
as  expected  from  the  theory.   Figure  2-6  shows  a  typical  case,  the 
Sepralyte  C-18  column,  with  acetonitrile  as  component  B  and  water  as 
component  A  of  the  mobile  phase.   The  exceptions  noted  are  the  same  as 
those  discussed  above.   The  free  energy  of  transfer  of  a  solute  per  unit 
area  can  be  calculated  by  utilizing  Figure  2-6  and  Table  2-5 

<Xsb  "  ZsA)k;T  =  -(3.689  x  10"2)  (592  cal/mol)  =  -21.8  cal/mol  A2   (2-8) 

and  is  approximately  constant  for  the  same  modifier  used.   The  above 
tests  confirm  the  predicted  dependence  of  solute  retention  on  cavity 
size,  and  the  chromatographic  method  provides  a  convenient  way  to 
experimentally  determine  solute-solvent  interaction  free  energies. 
According  to  equation  2-3,  the  slope  of  the  composition  plot 
divided  by  n  should  give  the  binary  interaction  constant  of  the  mobile 
phase  components,  XaB'  regardless  of  the  nature  of  the  solutes  and 
stationary  phases  used.   This  binary  interaction  constant  should  depend 
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y-mtercepts 
at  <J)B=1  from 
equation    2-3 
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Figure    2-5.    Plot    of   y-intercepts    at    (J>b=1    from  equation    2-3    versus 
van  der  Waals   volume,    Vw,    of   the   solutes    for   the 
Sepralyte   C-18    column   with    acetonitrile    as   modifier. 
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Table  2-4.      Regression  results  of  y-intercepts  at  <J>g=l  from  equation 
2-3  versus  the  van  der  Waals  volume,  Vw,  of  the  solutes 
for  all  the  columns . 


£T 


Column  Modifier         Slope (Ml Q21 

Sepralyte  C-2  Acetonitrile  -8.523  0.9542 

Sepralyte  C-4  Acetonitrile  -7.484  0.9344 

Sepralyte  C-8  Acetonitrile  -9.793  0.9112 

Sepralyte  C-18  Acetonitrile  -6.936  0.9178 

Sepralyte  C-4  Methanol  -12.73  0.9489 

Sepralyte  C-8  Methanol  -12.17  0.9136 

Hypersil  CDS  Acetonitrile  -377.2  0.7089 

Unisil  Q  C-18  Acetonitrile  -109.0  0.5035 

Ultrasphere  ODS  Methanol  -22.23  0.9993 

Ultrasphere  ODS  Ethanol  10.96  0.2696 

Ultrasphere  ODS  1-Propanol  95.35  0.1945 

Ultrasphere  ODS  Acetonitrile  -345.8  0.5705 

Silasorb  C-8  Methanol  721.2  0.2434 

Silasorb  C-8  Acetonitrile  629.9  0.2139 
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y-intercepts 
at  63=1  from 
eauation  2-3 
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Figure  2-6.  Plot  of  y-intercepts  at  <J>b=1  from  equation  2-3  versus 
hydrocarbonaceous  surface  area,  HSA,  of  the  solutes 
for  the  Sepralyte  C-18  column  with  acetonitrile  as 
modifier . 
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Table  2-5.      Regression  results  of  y-intercepts  at  Oa=l  from  equation 
2-3  vs  the  hydrocarbonaceous  surface  area,  HSA,  of  the 
solutes  for  all  the  columns. 


Column  Modifier         Slope  (>:102)  r.2 

Sepralyte  C-2  Acetonitrile  -4.383  0.9658 

Sepralyte  C-4  Acetonitrile  -3.957  0.9389 

Sepralyte  C-8  Acetonitrile  -5.118  0.9052 

Sepralyte  C-18  Acetonitrile  -3.689  0.9282 

Sepralyte  C-4  Methanol  -6.752  0.9460 

Sepralyte  C-8  Methanol  -6.378  0.9489 
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solely  on  the  nature  of  the  organic  modifier  in  the  mobile  phase  of 
RPLC.   By  comparing  the  values  of  the  slopes  of  Table  2-3,  they  are 
found  to  be  independent  of  the  solutes,  but  not  the  stationary  phase. 
There  is  a  factor  of  2  differentiating  the  slopes  from  a  short-chain 
stationary  phase  to  a  long-chain  stationary  phase,  such  as  Sepralyte  C-2 
to  the  Sepralyte  C-18  column.   This  variation  can  be  attributed  to  two 
factors.   First,  in  short-chain  stationary  phases,  surface  silanols  are 
more  accessible  to  solutes,  and  the  interactions  between  them  are  much 
larger  than  in  long-chain  phases.   Second,  an  alternate  retention 
mechanism,  adsorption  (Dill  1987) ,  should  dominate  in  short-chain 
stationary  phases.   Both  of  these  factors  are  neglected  in  the  original 
partition  model  (Dill  1987)  . 

Furthermore,  the  y-intercept  at  <|>B=1  divided  by  n  should  only  be 
dependent  on  the  nature  of  the  solute  and  organic  modifier.   The  data  of 
the  Sepralyte  columns  from  Table  2-4  confirm  this  prediction  since  the 
same  solutes  are  used  for  the  Sepralyte  columns.   The  large  differences 
of  the  values  in  column  3  of  Table  2-4  are  mainly  due  to  the  differences 
in  solutes  employed  with  different  stationary  phases. 

Johnson  and  coworkers  (Johnson  et  al.,  1986;  Dorsey  and  Johnson, 
1987)  have  shown  in  an  earlier  study  that  a  solvatochromic  dye  molecule, 
referred  to  as  ET(30),  can  be  used  to  probe  the  chemical  nature  of  the 
mobile  phase  and  its  strength  as  used  in  RPLC.   The  visible  spectral 
shift  of  the  dye  is  found  to  be  linearly  proportional  to  the  binary 
interaction  constant  %AB  of  the  mobile  phase  components  (Figure  2-7) . 
This  provides  a  justification  for  the  employment  of  ET(30)  as  a  probe 
for  the  study  of  the  mobile  phase,  since  it  appears  to  directly  measure 
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Figure  2-7.  Plot  slopes  from  Et(30)  plots  versus  slopes  from 
equation  2-3  for  the  Sepralyte  C-18  column  with 
acetonitrile  as  modifier. 
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Table  2-6.      Regression  results  of  slopes  from  equation  2-3  versus  the 
slopes  of  Erp(30)  plots  for  all  the  columns. 


r- 


Column  Modifier         Slope (xlQ2) 

Sepralyte  C-2  Acetonitrile  3.238  0.9996 

Sepralyte  C-4  Acetonitrile  2.699  0.9909 

Sepralyte  C-8  Acetonitrile  5.396  0.8796 

Sepralyte  C-18  Acetonitrile  4.047  0.9997 

Sepralyte  C-2  Methanol  14.67  0.9956 

Sepralyte  C-4  Methanol  9.315  0.9820 

Sepralyte  C-8  Methanol  11.03  0.9955 

Silasorb  C-8  Methanol  7.342  0.8369 

Silasorb  C-8  Acetonitrile  5.137  0.9987 

Unisil  Q  C-18  Acetonitrile  4.913  0.9303 

Hypersil  ODS  Acetonitrile  3.671  0.9187 

Ultrasphere  ODS  Methanol  8.055  0.9997 

Ultrasphere  ODS  Ethanol  6.166  0.8827 

Ultrasphere  ODS  1-Propanol  2.467  0.2630 

Ultrasphere  ODS  Acetonitrile  4.066  0.8702 
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the  free  energy  of  contact  between  components  A  and  B  of  the  mobile 
phase . 

For  the  Sepralyte  C-18  column  with  acetonitrile  as  modifier,  the 
slope  of  these  two  parameters  are  found  to  be 

(4.047  x  10"2)  x  (592  cal/mol)  =  23.96  units (cal/mol)       (2-9) 

Other  columns  are  observed  to  have  the  same  linearity  with  other  organic 
modifiers,  and  the  linear  regression  results  between  these  two 
parameters  are  listed  in  Table  2-6. 

Even  though  the  experimental  data  base  used  to  test  the 
theoretical  predictions  is  relatively  large,  the  range  of  mobile  phase 
compositions  is  still  limited  to  within  30  -  50%  variation  in  solvent 
composition.   The  linear  expression  proposed  by  Dill  (1987)  and  Dorsey 
and  Dill  (1989)  tested  in  this  work  is  found  to  hold  well  with  the  data 
base  used.   If  the  partition  model  and  other  assumptions  are  incorrect, 
nonlinearity  should  be  observed.   Since  good  linearity  is  found,  the 
slopes  and  y-intercepts  generated  by  plotting  the  composition  plot  of 
equation  2-3  should  reflect,  at  least  approximately,  the  binary 
interaction  constants  of  the  solutes  and  the  solvents,  and  can  be  used 
to  calculate  the  solute-solvent  interaction  free  energies  in  RPLC. 


CHAPTER  III 

CHARACTERIZATION  OF  THE  RETENTIVITY  OF 

REVERSED  PHASE  LIQUID  CHROMATOGRAPHY  COLUMNS 

Introduction 

There  is  an  immense  number  of  RPLC  columns  available  on  the  market 

today.   The  majority  of  them  are  made  from  Cg  or  C^g  functional  groups, 

but  the  use  of  columns  having  other  functionalities,  such  as  phenyl  and 
cyano,  has  been  on  the  rise.   Recent  studies  done  on  commercial  RPLC 
columns  have  found  that  they  all  showed  significant  differences  in 
absolute  retention  and  selectivity,  a,  for  the  same  solute  and  mobile 
phase,  even  when  they  are  all  labelled  as  C]o  columns  (Goldberg  1982; 

Sander  and  Wise,  1988)  .   This  variability  of  RPLC  columns  is  largely  due 
to  the  difference  in  the  starting  silica  as  well  as  the  bonding  reaction 
(Dorsey  and  Dill,  1989) .   Since  the  stationary  phase  in  RPLC  has  been 
shown  to  have  ample  effects  on  solute  retention  and  selectivity  (Sander 
1988;  Dill  1987;  Sentell  and  Dorsey,  1989a;  Sentell  and  Dorsey,  1989b), 
the  variability  has  caused  practical  chromatographers  many  difficulties 
in  choosing  the  best  column  to  develop  optimal  separations.   Developed 
methods  are  hard  to  transfer  unless  the  brand  and  manufacturer  of  the 
column  are  specified. 

Many  studies  have  been  performed  by  investigators  to  characterize 
the  retention  behavior  of  RPLC  columns.   Smith  (1982a  and  1982b) 
employed  a  homologous  series  of  alkylarylketones  to  develop  a  retention 
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index  with  a  set  of  reference  compounds  to  define  the  retention 
performance  of  RPLC  columns  made  from  different  manufacturers  and 
functionalities.   A  different  set  of  constants  for  each  column  is 
obtained  for  every  different  mobile  phase  composition  used.   These 
column  retention  constants  are  not  very  useful  since  new  calibration  on 
the  column  has  to  be  done  when  mobile  phase  composition  is  changed. 
Antle  and  coworkers  (Antle  and  Snyder,  1985;  Antle  et  al.,  1985)  used 
gradient  elution  theory  to  characterize  RPLC  columns  according  to  their 
solvophobic  retention.   They  examined  columns  produced  from  the  same 
base  silica  but  having  different  bonded  functionalities  such  as  C^g,  cq, 
phenyl,  Ct,  or  cyano  groups.   A  large  variety  of  test  solutes  having 
very  diverse  chemical  structure  were  used  in  their  study.   They  combined 
the  volume  phase  ratio,  <£>,    and  the  polarity  of  a  column  into  an 
effective  phase  ratio,  J',  which  gives  the  retentivity  of  a  column.   A 
reference  and  a  standard  column  were  used  to  acquire  the  relative  J' 
value  of  all  columns  studied.   Their  results  revealed  that  the  J'  value 
of  the  columns  are  in  the  order  of  C]^  >  C8  >  phenyl  >  C1    >  cyano.   They 
concluded  that  the  contribution  of  the  polarity  of  a  column  to  the 
retentivity  is  small  compared  to  the  phase  ratio  of  the  column.   Cooper 
and  Lin  (1986)  have  looked  at  the  retention  behavior  of  C3,  phenyl  and 
cyano  columns  using  three  solutes  chosen  from  Snyder's  selectivity 
triangle  (Snyder  1974;  1978) .   They  found  that  the  polarity  of  these 
columns  is  in  the  order  of  phenyl  >  cyano  >  Cg,  and  the  overall 
retentivity  of  these  columns  is  dominated  by  the  phase  ratio  of  the 
columns .   Other  researchers  have  used  chemometric  and  factor  analyses  to 
characterize  commercial  columns  (Delaney  et  al.,  1987;  Chretien  et  al., 
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1986;  Walczak  et  al.,  1987),  but  only  qualitative  results  were  obtained 
in  these  studies.   Walczak  et  al.  (1987)  concluded  that  carbon  loading, 
nature  of  the  organic  ligands,  and  the  accessibility  of  the  surface 
silanol  groups  are  the  main  factors  governing  the  retentivity  and 
selectivity  of  the  columns.   Delaney  et  al.  (1987)  summarized  that  their 
classification  of  RPLC  columns  from  chemometric  results  agreed  well  with 
a  qualitative  scheme  developed  by  a  liquid  chromatography  specialist. 

The  present  study  provides  a  simple  method  to  characterize  the 
retentivity  of  commercial  RPLC  columns.   Binary  mobile  phases  of  water 
and  an  organic  modifier  under  isocratic  conditions  were  used  throughout 
the  study  to  avoid  lengthy  equation  derivation,  and  to  keep  experimental 
parameters  simple.   No  reference  or  standard  column  was  needed,  and 
solutes  of  many  different  types  were  employed  so  that  the  results  of  the 
study  should  be  applicable  to  all  solutes. 

Experimental  Section 
All  the  retention  measurements  were  obtained  either  with  a 
Spectra-Physics  SP8700  ternary  proportioning  LC  system  (Spectra-Physics, 
San  Jose  CA,  USA)  or  with  two  Spectroflow  400  pumps  (ABI  Analytical, 
Kratos  Division,  Ramsey,  NJ,  USA) .   A  Valco  liquid  chromatography 
injection  valve  (Valco  Instrument  Company,  Houston,  TX,  USA)  and  a  20  |ll 
sample  loop  were  used  to  inject  the  solutes.  The  detection  system  was 
either  a  Spectroflow  783  absorbance  detector/gradient  controller  (ABI 
Analytical,  Kratos  Division,  Ramsey,  NJ,  USA)  operated  at  254  nm  or  a 
Spectroflow  757  absorbance  detector  (Kratos  Analytical  Instruments, 
Ramsey,  NJ,  USA)  also  operated  at  254  nm.   When  the  two  Spectroflow  400 
pumps  were  used,  a  Rainin  dynamax  dual  chamber  mixer  (Rainin  Instrument 
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Company,  Inc.,  Woburn,  MA,  USA)  was  placed  before  the  injection  valve 
for  better  mixing  of  the  binary  solvents.   All  the  Zorbax  columns  were 
commercially  available  from  Du  Pont  (E.  I.  du  Pont  de  Nemours  and 
Company,  Wilmington,  DE,  USA) ,  the  B  &  J  columns  from  Burdick  &  Jackson 
Laboratories,  Inc.  (Baxter  Healthcare  Corp.,  Burdick  &  Jackson  Division, 
Muskegon,  MI,  USA) ,  and  the  Ultrasphere  column  from  Beckman  (Beckman 
Instruments,  Inc.,  Fullerton,  CA,  USA).   The  high  density  column  was 
synthesized  and  packed  in  our  laboratories  (Sentell  1987) .   Some  of  the 
properties  of  these  columns  as  supplied  by  the  manufacturer  are  listed 
in  Table  3-1.   A  Fisher  Recordall  Series  5000  (Fisher  Scientific, 
Pittsburgh,  PA,  USA)  strip  chart  recorder  was  used  to  record  all  the 
chromatographic  peaks.   The  solvents  used  were  HPLC  grade  from  Fisher 
Scientific  (Fisher  Scientific,  Pittsburgh,  PA,  USA) ,  and  the  water  was 
filtered  through  a  Barnstead  NANOpure  II  system  (Barnstead  Company, 
Dubuque,  IA,  USA)  before  being  used.   The  columns  were  maintained  at 
30°C  by  a  water  jacket  and  a  Haake  Dl  circulator  (Haake,  Dieselstrasse, 
West  Germany) . 

A  total  of  26  solutes  was  used  and  they  are  listed  in  Table  3-2. 
They  were  from  various  suppliers  and  were  used  without  further 
purification.   The  solutes  were  grouped  into  mixtures  containing  no  more 
than  5  solutes  based  on  their  retention  data  when  they  were  injected 
individually  at  the  beginning  of  the  study.   Peak  assignments  were 
accomplished  by  identifying  the  peak  area  and  elution  order.   The  void 
volume  of  each  column  was  determined  by  the  injection  of  water  at  mobile 
phase  composition  of  75%  water  and  25%  organic  modifier.   The  flow  rate 
was  maintained  at  1  ml/min  except  for  the  Burdick  &  Jackson  C-18  column 
which  was  used  with  a  flow  rate  of  1.5  ml/min.   All  the  retention  data 
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Table  3-1. 


Properties  of  the  RPLC  columns  as  supplied  by  the 
manufacturers . 


Column 
Zorbax  ODS 

Zorbax  C-8 
Zorbax  phenyl 
Zorbax  TMS 
Zorbax  cyano 
B  &  J  C-8 
B  &  J  C-18 


Ultrasphere  C-8    Cg 


Ligand 
c18 

c8 

Phenylpropyl 

Cl 
Cyanopropyl 

C8 

c18 

Cc 


Surface  area(m2/g)    \   Carbon  loading 


340 
330 
300 
320 
320 
250 
250 
200 


20 

10 

5 

5 

5 

8 

13 

12 
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Table  3-2. 


List  of  test  solutes  used  in  this  study. 


Methyl  paraben 

Propyl  paraben 

Cortisone 

Toluene 

Benzyl  Alcohol 

Propachlor 

Tri-p-tolyl  phosphate 

Fluorobenzene 

Dimethyl  phthalate 

Butyl  benzyl  phthalate 

1-Methyl  phenanthrene 

p-Terphenyl 

1,1,4,4, -Teraphenyl-1, 3-butadiene 


Ethyl  paraben 
Butyl  paraben 
Corticosterone 
Hexyl fluorobenzene 
o-Nitrophenol 
1-Methyl  naphthalene 
Methyl  benzyl  amine 
Chloropropham 
Diethyl  phthalate 
Dibutyl  phthalate 
Dioctyl  phthalate 
Chrysene 
Octadecanophenone 
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were  obtained  from  averaging  at  least  two  measurements  of  a  solute.   The 
ET(30)  polarity  values  for  the  different  mobile  phase  compositions  were 
calculated  from  the  quadratic  relationship  between  percent  organic 
modifier  and  the  ET(30)  values  reported  by  Dorsey  and  Johnson  (1987). 

The  regression  calculations  were  all  done  by  using  the  StatWorks 
(Heydon  and  Son  Inc.,  Philadelphia,  PA,  USA)  program  on  a  Macintosh  SE 
(Apple  Computer,  Inc.,  Cupertino,  CA,  USA)  microcomputer. 

Results  and  Discussion 
Since  the  composition  of  the  mobile  phase  plays  a  primary  role  in 
the  retention  of  solutes  in  RPLC,  it  is  important  to  have  a  method  to 
account  for  the  solute  retention  contributed  by  the  mobile  phase  before 
the  retentivity  of  the  stationary  phase  can  be  explored.   Many  workers 
have  used  solvatochromic  measurements  to  measure  the  polarity  of  the 
mobile  phase  used  in  RPLC,  such  as  the  "Z"  scale  (Kosower  1958)  ,  the  71* 
multiparameter  scale  (Taft  and  Kamlet,  1976;  Kamlet  et  al.,  1983;  Taft 
et  al.,  1985;  Sadek  et  al.,  1985),  and  the  ET(30)  scale  (Johnson  et  al., 
1986;  Dorsey  and  Johnson,  1987).   Johnson  and  coworkers  (Johnson  et  al., 
1986;  Dorsey  and  Johnson,  1987)  have  shown  that  the  In  k '  of  a  solute  in 
RPLC  is  linearly  associated  with  the  E^(30)  polarity  scale  and  has 

better  linearity  than  the  volume  percent  of  organic  modifier.   Their 
relationship  can  be  expressed  as 

In  k'  =  m  (  ET (30)  )  +  c  (3-1) 
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where  Et>(30)  is  the  polarity  value  in  kcal/mol  of  the  mobile  phase 
measured  by  the  ET(30)  scale,  and  m  and  c  are  the  slope  and  y-intercept 
of  the  linear  regression.   The  ET(30)  values  are  found  to  have  an 

excellent  quadratic  relationship  with  most  common  binary  mobile  phases 
used  in  RPLC  (Dorsey  and  Johnson  1987) .   Moreover,  the  slope  m  has  been 
shown  in  Chapter  II  to  be  directly  proportional  to  the  size  of  the 
solutes  and  the  binary  interaction  constant  of  the  mobile  phase.   We 
selected  to  employ  the  slope  m  from  equation  3-1  as  the  descriptor  for 
the  mobile  phase  contribution  to  solute  retention. 
Molecular  Descriptor  Approach 

Although  the  mobile  phase  effect  can  be  accounted  for  by  using  the 
slope  in  equation  3-1,  the  cavity  created  in  the  stationary  phase  also 
is  a  major  contribution  to  solute  retention  according  to  the  partition 
model.   The  first  approach  we  took  in  this  study  was  to  search  for  a 
molecular  descriptor  of  the  test  solutes  that  can  delineate  the  cavity 
in  the  stationary  phase,  and  use  it  to  plot  against  slope  m.   The 
resulting  slope  of  this  plot  should  only  be  a  function  of  the 
retentivity  of  the  column.   We  first  predicted  that  a  linear 
relationship  should  occur  between  the  slopes  from  equation  3-1  and  a 
molecular  size  descriptor.   It  has  been  shown  in  the  literature  that  a 
linear  dependence  exists  between  In  k'  of  a  solute  and  some  molecular 
size  descriptors  of  the  solute  such  as  the  van  der  Waals  volume,  Vw, 
(Hanai  and  Hubert,  1984a;  Smith  1981;  Jinno  and  Kawasaki,  1983a)  and  the 
molecular  connectivity  index,  %l     (Wells  et  al.,  1982;  1986).   Slopes 
from  equation  3-1  were  plotted  against  both  the  van  der  Waals  volume  and 
the  molecular  connectivity  index  (Figures  3-1  and  3-2)  and  the 
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Slopes  from 
equation  3-1 
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Figure  3-1.  Plot  of  slopes  from  equation  3-1  against  the  van  der 
Waals  volume,  Vw,  of  the  test  solutes  for  the  Zorbax 
TMS  column  using  acetonitrile  as  modifier. 
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Table  3-3, 


Regression  results  of  graphs  of  slopes  from  equation  3-1 
against  van  der  Waals  volume,  Vw,  of  the  test  solutes. 


Column 


Modifier      SlopexdO3) 


y-inter. 


Zorbax  ODS  ACN 

Zorbax  C-8  ACN 

Zorbax  phenyl  ACN 

Zorbax  TMS  ACN 

Zorbax  cyano  ACN 


4.08 
3.82 
4.27 
9.25 
6.98 


0.845  0.533 

1.03  0.534 

0.593  0.484 

0.336  0.861 

0.310  0.472 
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Slopes  from 
equation  3-1 


connectivity  index 


Figure  3-2.  Plot  of  slopes  from  equation  3-1  against  the  molecular 
connectivity  index  of  the  test  solutes  for  the  Zorbax 
TMS  column  using  acetonitrile  as  modifier. 
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Table  3-4. 


Regression  results  of  graphs  of  slopes  from  equation  3-1 
against  molecular  connectivity  index  of  the  test  solutes, 


Column 


Modifier      SlopexdO2)        y-inter. 


Zorbax  ODS  ACN 

Zorbax  C-8  ACN 

Zorbax  phenyl  ACN 

Zorbax  TMS  ACN 

Zorbax  cyano  ACN 


6.21 
6.54 
7.14 
16.7 
11.5 


0.992  0.401 

1.14  0.417 

0.723  0.379 

0.563  0.829 

0.531  0.354 
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regression  results  of  these  graphs  on  five  Zorbax  columns  are  shown  in 
Table  3-3  and  3-4  respectively. 

The  coefficients  of  determination,  r  ,  of  these  regressions  are 
all  well  below  0 . 9  meaning  that  there  is  no  significant  correlation 
between  the  slopes  from  equation  3-1  and  the  two  molecular  size 
descriptors  of  the  solutes  that  we  have  chosen.   One  of  the  reasons  for 
the  breakdown  of  our  initial  hypothesis  is  probably  because  the 
molecular  size  descriptors  do  not  account  for  all  the  interactions  among 
the  solutes  and  the  stationary  phase.   Despite  Chapter  II  showing  that 
the  cavity  opened  in  the  stationary  phase  is  proportional  to  the  size  of 
the  solute,  there  are  other  chemical  interactions  such  as  polarity, 
dipole  moment,  and  hydrogen  bonding  ability  between  the  solute  and  the 
stationary  phase.   These  interactions  cannot  be  completely  summed  by  one 
single  molecular  size  descriptor. 
In  kw  Approach 

The  next  approach  we  took  was  to  seek  a  parameter  that  will 
accurately  include  all  the  interactions  between  the  solute  and  the 
stationary  phase.   The  logarithmic  capacity  factor  of  a  solute  at  100% 
water,  In  kw,  has  been  shown  to  have  great  correlation  to  the  logarithm 

of  the  water/octanol  partition  coefficient,  log  PQ/W,    of  a  solute 
(Miyake  et  al.,  1988;  Braumann  1986;  Braumann  et  al.,  1987;  Minick  et 
al.,  1987)  .   The  linear  relationship  shows  that  both  In  kw  and  log  P0/w 

are  measuring  a  similar  partition  process.   This  suggests  that  In  kw 
measures  the  driving  force  of  the  partition  of  a  solute  into  the 
stationary  phase,  and  is  the  parameter  that  measures  all  the 
interactions  between  the  solute  and  the  stationary  phase. 
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The  In  kw  of  a  solute  is  difficult  to  obtain  experimentally  in 
RPLC  due  to  long  retention  times  and  poor  peak  shapes  arising  from  slow 
mass  transfer  with  pure  water  as  the  mobile  phase.   There  is  a  popular 
belief  that  the  logarithm  of  the  capacity  factor,  In  k',  of  a  solute  is 
linearly  related  to  the  percent  by  volume  of  organic  modifier  in  the 
mobile  phase  (Snyder  et  al.,  1979;  Dolan  et  al.,  1979;  Schoemakers  et 
al.,  1979;  Antle  et  al.,  1985;  Baty  and  Sharp,  1988).    Hence,  In  kw  is 
often  estimated  by  extrapolation  from  a  linear  regression  of  In  k' 
versus  percent  organic  modifier  back  to  0%  organic  modifier  (Reymond  et 
al.,  1987;  Braumann  1986;  Braumann  et  al.,  1987;  Baty  and  Sharp,  1988). 

In  kw  =  S  (%organic)  +  c'  (3-2) 

Unfortunately,  this  kind  of  extrapolation  has  been  found  to  be  not  as 
reliable  as  using  the  ET(30)  polarity  scale  (Michels  1989).   The  ET(30) 
polarity  value  of  pure  water  is  63.11  kcal/mol,  and  therefore  In  kw  of  a 
solute  can  be  approximated  by  substituting  63.11  kcal/mol  into  equation 
3-1. 

Since  the  slopes  from  equation  3-1  and  the  In  kw  of  the  solutes 

summarize  the  retention  effects  due  to  the  mobile  phase  and  the  solute 
interactions  with  the  stationary  phase  respectively,  the  slope  of  a  plot 
of  these  two  parameters  for  a  column  should  be  independent  of  the  two 
effects,  and  unveil  the  retentivity  of  the  column  (Figure  3-3) .   The 
regression  results  of  these  plots  for  the  columns  studied  are  presented 
in  Table  3-5.   The  95%  confidence  interval  of  the  slopes  of  these 
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Figure  3-3.  Plot  of  slopes  from  equation  3-1  against  the  In  kw  of 
the  test  solutes  for  the  Zorbax  TMS  column  using 
acetonitrile  as  modifier. 


Table  3-5.      Regression  results  of  slopes  from  equation  3-1  vs  In  kv 
of  the  test  solutes  for  all  the  columns. 


Column 

Modifier 

Zorbax 

ODS 

ACN 

0 

0835 

Zorbax 

C-8 

ACN 

0 

069 

Zorbax 

Phenyl 

ACN 

0 

104 

Zorbax 

TMS 

ACN 

0 

110 

Zorbax 

Cyano 

ACN 

0 

126 

Zorbax 

ODS 

MeOH 

0 

120 

Zorbax 

C-8 

MeOH 

0 

124 

Zorbax 

Phenyl 

MeOH 

0 

129 

Zorbax 

TMS 

MeOH 

0 

134 

Zorbax 

Cyano 

MeOH 

0 

130 

B  &  J  C-8 

MeOH 

0 

134 

B  &  J  C-18 

MeOH 

0 

114 

Ultrasphere  C-8 

MeOH 

0 

122 

High  Density 

MeOH 

0 

115 

y-inter.  £^ 

0.525  0.894 

0.727  0.834 

0.293  0.98 

0.349  0.982 

0.241  0.99 

0.320  0.977 

0.310  0.987 

0.271  0.993 

0.346  0.979 

0.362  0.96 

0.209  0.987 

0.326  0.945 

0.307  0.99 

0.314  0.923 


CJ^o.os' 

±0.004 

±0.035 

±0.005 

±0.007 

±0.005 

±0.003 

±0.004 

±0.003 

±0.006 

±0.006 

±0.009 

±0.012 

±0.004 

±0.013 


a)  95%  confidence  interval  of  the  slopes  of  the  regressions 
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regressions  are  calculated  (Sharaf  et  al . ,  1986;  Anderson  1987),  and 
they  are  listed  in  the  last  column  of  Table  3-5.   The  confidence 
interval  shows  that  almost  all  the  slopes  of  the  regressions  are 
statistically  different  from  one  another.   Good  linearity  appears 
between  these  two  parameters  with  the  r2  of  all  the  columns  tested  above 
0.9,  except  in  the  case  of  the  Zorbax  ODS  and  C-8  columns  with 
acetonitrile  as  modifier.   The  high  r2  confirms  that  the  slopes  from 
equation  3-1  and  In  kw  of  the  test  solutes  are  highly  correlated.   This 
type  of  slope  versus  In  kw  plot  is  not  new.   Several  groups  (Braumann  et 
al.,  1983;  Hammer  et  al.,  1982)  have  shown  that  plots  of  the  slope,  S, 
and  In  kw  from  equation  3-2  are  linear.   They  found  the  slopes  and 
y-intercepts  of  these  plots  formed  two  empirical  parameters  that  can  be 
employed  to  classify  solutes  into  different  groups.   These  solute  groups 
can  be  used  to  recommend  the  use  of  In  kw  for  the  estimation  of  log 
PQ/W.   Schoenmakers  et  al.  (1979)  found  good  linearity  when  methanol  is 
used  as  mobile  phase  modifier  but  not  with  acetonitrile  or  THF .   They 
have  shown  that  when  linearity  is  observed,  the  slopes  of  these  plots 
can  be  used  to  determine  the  shape  of  the  gradient  program.   Moreover, 
these  slopes  can  be  employed  to  predict  isocratic  capacity  factors  from 
a  simple  gradient  analysis  (Schoenmakers  et  al.,  1981).   Baty  and  Sharp 
(1988)  observed  good  correlation  between  S  and  In  kw  for  methanol, 
acetonitrile  and  THF  using  structurally  similar  solutes.   They  tried  to 
use  the  slopes  and  y-intercepts  of  these  plots  to  predict  capacity 
factors  of  the  solutes,  but  they  found  that  the  slopes  and  y-intercepts 
of  these  plots  vary  with  the  nature  of  the  organic  modifiers  and 
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columns.   Therefore,  the  chromatographic  system  and  the  solute  group 
have  to  be  defined  before  the  capacity  factors  can  be  predicted. 

The  general  retentivity  of  the  columns  in  this  study  for  a  given 
organic  modifier  is  found  to  be  inversely  proportional  to  the  slope  of 
the  slopes  versus  In  kw  plot.   The  usefulness  of  this  relative 
retentivity  scale  can  be  seen  by  plotting  the  linear  regressions  of  the 
five  Zorbax  columns  using  acetonitrile  as  modifier  on  one  graph  (Figure 
3-4) .   As  shown  in  Chapter  II,  the  slopes  from  equation  3-1  are  found  to 
be  proportional  to  the  size  of  the  cavity  created  by  the  solutes; 
therefore,  at  a  given  solute  cavity  size,  the  linear  regression  having 
the  smallest  slope  will  have  the  largest  In  kw  and  hence,  the  largest 
retentivity.   From  the  data  shown  in  Figure  3-4,  the  Zorbax  cyano  column 
has  the  largest  slope  followed  by  the  TMS  column,  phenyl  column,  ODS 
column  and  C-8  column;  retentivity  of  these  column  is  in  the  order  of 
C-8  >  ODS  >  phenyl  >  TMS  >  cyano.   The  Zorbax  C-8  column  seems  to  have  a 
higher  retentivity  than  the  Zorbax  ODS  column  when  acetonitrile  is  used 
as  modifier.   A  careful  examination  of  the  confidence  interval  of  these 
two  slopes  demonstrates  that  they  actually  have  no  statistical 
difference,  and  the  apparent  higher  retentivity  of  the  Zorbax  C-8  column 
could  very  well  be  an  experimental  artifact.   The  apparent  stronger 
retentivity  of  the  Zorbax  cyano  column  over  the  TMS  column  with  methanol 
as  modifier  can  be  attributed  to  the  hydrogen  bonding  between  the 
surface  silanols  and  the  methanol.   Cyano  columns  have  been  shown  to 
have  more  accessible  silanols  than  other  columns  (Cooper  and  Lin,  1986; 
Smith  and  Miller,  1989) .   Since  methanol  can  form  hydrogen  bonds  with 
these  silanols,  the  surface  of  cyano  columns  is  significantly  modified 
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Slopes    from 
equation    3-1 


10  20 

In  kw  of  solutes 


30 


Figure  3-4.  Plot  of  slopes  from  equation  3-1  against  the  In  kw  of 
the  test  solutes  for  all  the  Zorbax  columns  using 
acetonitrile  as  modifier. 
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by  the  methanol,  and  hence  its  apparent  retentivity  is  stronger  than  the 
TMS  column.   When  acetonitrile  is  used  as  modifier,  due  to  its  lack  of 
hydrogen  bonding  ability,  the  retentivity  of  the  TMS  column  is  found  to 
be  stronger  than  the  cyano  column.   Other  commercial  columns  can  be 
tested  using  the  same  procedures  and  compared  with  the  retentivity  of 
the  Zorbax  columns.   With  methanol  as  modifier,  the  B  &  J  C-18  column  is 
found  to  have  retentivity  larger  than  any  of  the  Zorbax  columns,  and  the 
B  &  J  C-8  column  is  found  to  have  retentivity  approximately  the  same  as 
the  Zorbax  TMS  column,  while  the  Ultrasphere  C-8  column  is  shown  to  have 
retentivity  between  the  Zorbax  ODS  and  C-8  columns. 

One  interesting  point  raised  by  many  researchers  is  the  importance 
of  phase  ratio  of  the  column  on  the  total  retentivity  of  the  column. 
The  significance  of  phase  ratio  on  retentivity  was  also  investigated  in 
this  study.   Phase  ratio  in  RPLC  can  be  calculated  using  the  volume 
ratio  of  the  stationary  phase  to  the  mobile  phase.   Although  much  work 
has  been  done  on  measuring  the  volume  of  the  mobile  phase,  Vm,  (Melander 
et  al.,  1983;  Smith  et  al.,  1986;  Engelhardt  et  al . ,  1984;  Hennion  and 
Rosset,  1988)  measurement  of  the  volume  of  the  stationary  phase,  Vs,  is 
more  ambiguous  (Jandera  et  al.,  1982;  Melander  et  al.,  1980;  Sander  and 
Field,  1980;  Slaats  et  al.,  1981;  McCormick  and  Karger,  1980a; 
Berendsen  et  al.,  1980) .   By  employing  the  Vs  and  Vm  calculation  method 
presented  by  Sentell  (1987),  the  phase  ratios  of  the  columns  studied 
were  obtained  and  are  listed  in  Table  3-6.   In  the  case  of  the  phenyl 
and  cyano  columns,  the  density  of  the  ligands  are  not  available  and  the 
Vs  of  the  columns  cannot  be  calculated.   The  phase  ratio  of  these 
columns  are  therefore  not  reported  in  Table  3-6.   With  the  exception  of 
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Table  3-6.      Phase  ratio  of  the  columns  in  this  study  calculated  using 
the  method  presented  by  Sentell  (1988)  . 


Column  Phase  ratio 

Zorbax  ODS  0.3  67 

Zorbax  C-8  0.226 

Zorbax  phenyl  unknown 

Zorbax  TMS  0.172 

Zorbax  cyano  unknown 

B  &  J  C-8  0.0619 

B  &  J  C-18  0.110 


Ultrasphere  C-8  0.380 

High  density  0.465 
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the  two  B  &  J  columns,  the  phase  ratio  of  the  columns  displays  the  same 
order  as  the  relative  retentivity  scale  of  the  columns.   The  high 
density  column,  which  has  the  largest  phase  ratio,  is  found  to  be  close 
to  having  the  highest  retentivity.   This  result  points  out  that  the 
phase  ratio  of  a  column  indeed  has  a  vital  contribution  to  the  overall 
retentivity  of  the  column. 

This  work  exhibits  a  simple  method  that  can  be  used  to  classify 
commercial  RPLC  columns  according  to  their  overall  retentivity.   The 
isocratic  approach  used  has  kept  the  data  interpretation  and 
manipulation  simple.   The  use  of  In  kw  in  the  procedure  is  justified  and 
is  found  to  be  very  appropriate.   The  phase  ratio  of  a  column  is  shown 
to  dominate  the  total  retentivity  of  the  column.   Since  such  a  large 
variety  of  test  solutes  were  employed  in  this  work,  the  relative 
retentivity  scale  based  on  this  study  should  be  applicable  to  almost  all 
solutes . 


CHAPTER  IV 
CONCLUSIONS  AND  FUTURE  WORK 

Conclusions 

The  research  that  has  been  discussed  in  the  previous  chapters  was 
performed  to  confirm  and  apply  the  partition  retention  mechanism  of 
RPLC.   The  recognition  of  the  partition  mechanism  as  the  correct 
retention  model  in  RPLC  is  very  important  because  it  can  help 
chromatographers  to  understand  the  processes  happening  in  RPLC.   Also, 
with  a  sound  fundamental  understanding  of  solute  retention,  more  robust 
separation  methods  can  be  developed. 

In  Chapter  II,  we  have  tested  the  partition  model  with  a  large 
data  base.   We  found  that  the  relationships  between  solute  capacity 
factors  and  mobile  phase  compositions  are  in  good  agreement  with  the 
prediction  from  the  partition  model  (Dill  1987;  Dorsey  and  Dill,  1989) 
with  only  very  minor  discrepancies.   The  partition  model  proposed  by 
Dill  (1987)  is  based  on  simple  lattice  models  and  the  interphase  model 
of  the  stationary  phase,  but  the  results  of  our  tests  are  still  very 
good.   With  a  more  explicit  treatment  of  the  simple  assumptions,  the 
theory  may  be  able  to  account  for  the  small  discrepencies .   The 
molecular  origins  of  these  composition  plots  are  also  clarified  in  our 
tests.   The  slopes  of  these  plots  are  found  to  be  a  function  of  the 
binary  interaction  parameter  of  the  mobile  phase  components,  and  the 
y-intercepts  at  $^=1   are  proportional  to  the  free  energy  of  transfer  of 
the  solutes.   The  solute  size  plots  provide  a  method  to  calculate  the 
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solute-solvent  interaction  free  energies,  and  an  explanation  to  the 
widely  observed  relationship  between  capacity  factors  and  solute  sizes 
in  RPLC  (Hanai  and  Hubert,  1984b;  Jinno  and  Kawasaki,  1983a;  1983b;  Feng 
et  al.,  1988;  Mockel  et  al.,  1987).   Moreover,  the  slopes  and 
y-intercepts  at  4>b=1  of  the  composition  plots  predicted  by  the  theory 
give  a  simple  and  reliable  method  to  estimate  solute-solvent  interaction 
free  energies.   These  free  energies  are  difficult  to  obtain 
experimentally  and  are  seldom  reported  in  the  literature.   With  this  new 
approach,  these  interaction  free  energies  should  be  much  easier  to 
access  due  to  the  availability  of  chromatographic  data. 

The  partition  model  is  put  to  practice  in  Chapter  III.   We 
employed  the  partiton  model  as  our  basis  to  develop  a  relative 
retentivity  scale  for  RPLC  columns.   In  the  partiton  model,  the  entire 
stationary  phase  contribution  to  retention  is  summed  in  the  In  kw  term 
of  the  solutes.   Although  this  assumption  is  an  oversimplification  of 
the  stationary  phase  effect,  it  has  been  shown  that  In  kw  is  a  very 

effective  parameter  to  estimate  the  solute  partition  coefficient  (Miyake 
et  al.,  1988;  Braumann  1986;  Braumann  et  al.,  1987;  Minick  et  al., 

1989) .   We  used  the  In  kw  of  the  solutes  and  the  slopes  from  the  In  k' 
versus  ET(30)  values  plots  to  obtain  the  relative  retentivity  values  of 

the  RPLC  columns.   The  correlation  between  these  two  parameters  is  shown 
to  be  good,  and  this  further  supports  the  partition  model  as  the 
dominant  solute  retention  mechanism.   The  retentivity  scale  is  found  to 
be  in  good  agreement  with  the  literature  (Antle  and  Snyder,  1985;  Antle 
et  al . ,  1985).   Also,  the  experimental  procedures  are  kept  relatively 
simple,  and  therefore  this  retentivity  scale  should  be  useful  in 
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classifying  the  retentivity  of  commercial  RPLC  columns,  and  should  help 
practical  chromatographers  to  select  the  best  column  for  their 
applications.   Since  test  solutes  of  various  chemical  structures  are 
employed  in  our  study,  the  retentivity  scale  should  be  applicable  to  all 
solutes.   The  phase  ratio  of  the  columns  is  found  to  play  a  major  role 
in  determining  the  retentivity  of  the  column. 

Future  Work 

Although  the  data  base  that  is  used  in  Chapter  II  to  test  the 
predictions  from  the  partition  model  is  rather  large,  it  lacks  a  wide 
range  in  mobile  phase  compositions  especially  at  the  extreme 
compositions.   Most  of  the  data  sets  have  mobile  phase  composition  range 
between  30  -  50%  organic  modifier.   A  test  of  the  predictions  at  extreme 
mobile  phase  compositions  such  as  mobile  phases  having  a  few  percent  of 
organic  modifier  would  be  interesting,  since  it  will  unveil  how  well  the 
theory  holds  under  extreme  mobile  phase  composition.   The  stationary 
phase  has  been  reported  to  take  on  a  different  structure  as  the 
composition  of  the  mobile  phase  changes  (Yonker  et  al.,  1982a;  1982b; 
McCormick  and  Karger,  1980a;  1980b;  Carr  and  Harris,  1986;  1987;  McNally 
and  Rogers,  1985) .   The  changing  in  the  structure  of  the  stationary 
phase  implies  that  the  retention  mechanism  may  switch  from  partition  to 
other  forces.   By  testing  the  theory  at  extreme  compositions,  the  mobile 
phase  composition  that  can  significantly  change  the  stationary  phase 
structure  may  be  observed  as  the  relationship  becomes  nonlinear. 

Another  test  that  is  useful  in  testing  the  validity  of  the 
partition  model  is  to  use  retention  data  from  stationary  phases  having 
functionalities  other  than  alkyl  chains,  such  as  phenyl  and  cyano 
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groups.   The  data  from  these  stationary  phases  should  give  information 
on  the  retention  mechanism  in  these  different  phases.   If  linearity 
exists  in  these  columns,  partition  should  be  the  dominate  force.   If  a 
nonlinear  relationship  occurs,  a  retention  mechanism  other  than 
partition  should  be  considered  for  these  columns. 

The  experimental  procedures  for  obtaining  the  relative 
retentiviity  scale  in  Chapter  III  is  simple;  nevertheless  it  is  rather 
time  consuming  since  isocratic  data  of  more  than  20  solutes  have  to  be 
collected.   A  possible  time  saving  step  is  to  find  a  few  representative 
solutes  so  as  to  eliminate  the  useage  of  all  2  6  test  solues .   The  Snyder 
selectivity  triangle  (Snyder  1974;  1978)  provides  a  good  starting  place 
to  look  for  representative  solutes. 

Moreover,  the  present  study  only  involved  monomeric  coverage 
columns.   In  order  to  expand  the  practical  utility  of  the  retentivity 
scale,  stationary  phases  other  than  monomeric  coverage  should  also  be 
tested  to  find  out  if  they  fit  in  the  retentivity  scales.   Since 
polymeric  phases  have  been  shown  to  have  better  selectivity  for  PAH's 
than  monomeric  phases  (Sander  and  Wise,  1984),  and  more  and  more 
polymeric  phases  are  available,  there  is  a  need  to  classify  and  compare 
them  with  the  monomeric  phases.   By  using  polymeric  columns,  the 
applicability  of  our  retentivity  scale  can  be  tested.   Also,  since  the 
retentivity  scale  is  based  on  the  partition  model,  a  failure  in 
classifying  polymeric  columns  with  the  retentivity  scale  may  suggest 
that  a  retention  mechanism  other  than  partition  is  the  dominate  factor 
in  polymeric  phases.   The  surface  coverage  of  polymeric  phases  is 
usually  higher  than  monomeric  phases.   It  has  been  pointed  out  that 
these  reported  values  may  not  give  the  true  surface  coverage  for 
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polymeric  phases  (Sander  and  Wise,  1984;  Dorsey  and  Dill,  1989),  because 
during  the  polymerization  process,  it  is  very  likely  that  a  reactive 
silane  molecule  will  bond  at  a  point  away  from  the  silica  surface.   The 
validity  of  these  reported  values  may  be  revealed  by  combining  the  phase 
ratio  and  the  retentivity  values  of  these  polymeric  columns. 


APPENDIX  A 
CHROMATOGRAPHIC  RETENTION  DATA 
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APPENDIX  B 
SLOPES  FROM  EQUATION  3-1  AND  In  kv 


Column:  Zorbax  ODS 
Mobile  phase:  ACN/water 
Reference:  This  work 

■Solute 

Methyl  Paraben 

Ethyl  Paraben 

Propyl  Paraben 

Butyl  Paraben 

Cortisone 

Corticosterone 

Toluene 

Hexylf luorobenzene 

Benzyl  Alcohol 

o-Nitrophenol 

Propachlor 

1-Methyl  naphthalene 

Tri-p-toyly  phosphate 

Methyl  benzyl  amine 

Fluorobenzeze 

Chloropropham 

Dimethyl  phthalate 

Diethyl  phthalate 

Dibutyl  phthalate 

1-Methyl  phenanthrene 

Butyl  benzyl  phthalate 

Dioctyl  phthalate 

p-Terphenyl 

Chrysene 


in_Xv 


0.95 

5.14 

1.21 

7.56 

1.25 

8.09 

1.39 

9.52 

1.39 

7.82 

1.28 

8.08 

1.06 

8.50 

1.39 

10.22 

0.65 

3.25 

0.98 

6.66 

1.19 

8.40 

1.38 

11.39 

2.00 

16.22 

1.03 

6.30 

1.20 

8.63 

1.49 

10.93 

1.41 

8.99 

1.41 

9.89 

1.61 

12.91 

1.50 

12.77 

1.70 

13.69 

1.89 

18.28 

1.66 

14.76 

1.49 

13.23 
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Column:  Zorbax  C-8 
Mobile  phase:  ACN/water 
Reference:  This  work 

Solute  Slope  from  equation  3-1 


Methyl  Paraben 

Ethyl  Paraben 

Propyl  Paraben 

Butyl  Paraben 

Cortisone 

Corticosterone 

Toluene 

Hexylf luorobenzene 

Benzyl  Alcohol 

o-Nitrophenol 

Propachlor 

1-Methyl  naphthalene 

Tri-p-toyly  phosphate 

Methyl  benzyl  amine 

Fluorobenzeze 

Chloropropham 

Dimethyl  phthalate 

Diethyl  phthalate 

Dibutyl  phthalate 

1-Methyl  phenanthrene 

Butyl  benzyl  phthalate 

Dioctyl  phthalate 

p-Terphenyl 

Chrysene 


1.28 

7.18 

1.24 

7.46 

1.57 

10.18 

1.56 

10.55 

1.56 

8.05 

1.54 

9.29 

1.27 

9.15 

1.36 

10.43 

0.82 

3.84 

1.31 

8.77 

1.34 

9.07 

1.39 

11.22 

1.72 

14.45 

1.27 

8.25 

1.25 

8.69 

1.50 

11.07 

1.40 

8.75 

1.38 

9.79 

1.60 

13.38 

1.43 

12.15 

1.68 

13.42 

2.20 

20.84 

1.67 

13.89 

1.56 

13.25 
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Column:  Zorbax  phenyl 
Mobile  phase:  ACN/water 
Reference:  This  work 

Solute 


Methyl  Paraben 

Ethyl  Paraben 

Propyl  Paraben 

Butyl  Paraben 

Cortisone 

Corticosterone 

Toluene 

Hexylf luorobenzene 

Benzyl  Alcohol 

o-Nitrophenol 

Propachlor 

1-Methyl  naphthalene 

Tri-p-toyly  phosphate 

Methyl  benzyl  amine 

Fluorobenzeze 

Chloropropham 

Dimethyl  phthalate 

Diethyl  phthalate 

Dibutyl  phthalate 

1-Methyl  phenanthrene 

Butyl  benzyl  phthalate 

Dioctyl  phthalate 

p-Terphenyl 

Chrysene 

1,1,4,4, -Tetraphenyl- 

1, 3-butadiene 

Octadecanophenone 


In  k^, 


0.81 

4.68 

0.92 

5.60 

1.07 

7.23 

1.20 

7.93 

1.00 

5.72 

0.99 

6.36 

0.92 

6.66 

1.04 

7.63 

0.56 

2.87 

0.81 

5.33 

0.98 

6.80 

1.20 

9.33 

1.60 

12.98 

0.79 

5.02 

0.90 

6.33 

1.24 

8.76 

0.87 

5.65 

1.01 

7.14 

1.44 

11.08 

1.40 

10.85 

1.77 

13.00 

2.11 

17.16 

1.81 

14.23 

1.67 

12.89 

2.12 

17.79 

2.47 

20.88 
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Column:  Zorbax  TMS 
Mobile  phase:  ACN/water 
Reference:  This  work 

Solute 

Methyl  Paraben 

Ethyl  Paraben 

Propyl  Paraben 

Butyl  Paraben 

Cortisone 

Corticosterone 

Toluene 

Hexylf luorobenzene 

Benzyl  Alcohol 

o-Nitrophenol 

Propachlor 

1-Methyl  naphthalene 

Tri-p-toyly  phosphate 

Methyl  benzyl  amine 

Fluorobenzeze 

Chloropropham 

Dimethyl  phthalate 

Diethyl  phthalate 

Dibutyl  phthalate 

1-Methyl  phenanthrene 

Butyl  benzyl  phthalate 

Dioctyl  phthalate 

p-Terphenyl 

Chrysene 

1,1,4,4, -Tetraphenyl- 

1, 3-butadiene 

Octadecanophenone 


In  kv 


1.01 

5.24 

1.02 

5.97 

1.17 

6.94 

1.28 

8.18 

1.18 

6.07 

1.19 

6.90 

1.03 

6.60 

1.12 

7.58 

0.51 

1.89 

0.59 

3.68 

1.21 

7.56 

1.31 

9.17 

2.20 

15.84 

0.77 

4.25 

0.99 

6.11 

1.34 

8.97 

0.90 

5.24 

1.06 

7.00 

1.71 

12.52 

1.60 

11.38 

1.81 

13.23 

2.39 

18.83 

1.95 

14.06 

1.72 

12.55 

2.27 

17.26 

2.52 

21.04 
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Column:  Zorbax  cyano 
Mobile  phase:  ACN/water 
Reference:  This  work 

Solute  Slope  from  equation  3-1      In  k 

Methyl  Paraben 

Ethyl  Paraben 

Propyl  Paraben 

Butyl  Paraben 

Cortisone 

Corticosterone 

Toluene 

Hexylf luorobenzene 

Benzyl  Alcohol 

o-Nitrophenol 

Propachlor 

1-Methyl  naphthalene 

Tri-p-royly  phosphate 

Methyl  benzyl  amine 

Fluorobenzeze 

Chlorcpropham 

Dimethyl  phthalate 

Diethyl  phthalate 

Dibutyl  phthalate 

1-Methyl  phenanthrene 

Butyl  benzyl  phthalate 

Dioctyl  phthalate 

p-Terphenyl 

Chrysene 

1,1,4,4, -Tetraphenyl- 

1, 3-butadiene 

Octadecanophenone 


0.67 

3.18 

0.72 

3.88 

1.04 

5.83 

1.14 

6.95 

0.93 

4.39 

0.95 

5.08 

0.92 

5.73 

1.12 

7.18 

0.38 

1.19 

0.57 

3.40 

0.94 

5.51 

1.21 

8.26 

2.11 

14.16 

0.52 

2.69 

0.72 

4.33 

1.13 

7.21 

0.99 

5.28 

1.10 

6.72 

2.04 

13.74 

1.73 

11.68 

1.60 

11.28 

2.55 

17.93 

1.90 

13.91 

1.75 

12.54 

2.34 

16.68 

2.69 

19.77 
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Column:  Zorbax  ODS 
Mobile  phase:  MeOH/water 
Reference:  This  work 

golute  Slope  from  equation  3-1      In  Kw 

Methyl  Paraben 

Ethyl  Paraben 

Propyl  Paraben 

Butyl  Paraben 

Cortisone 

Corticosterone 

Toluene 

Hexylf luorobenzene 

Benzyl  Alcohol 

o-Nitrophenol 

Propachlor 

1-Methyl  naphthalene 

Tri-p-toyly  phosphate 

Methyl  benzyl  amine 

Fluorobenzeze 

Chloropropham 

Dimethyl  phthalate 

Diethyl  phthalate 

Dibutyl  phthalate 

1-Methyl  phenanthrene 

Butyl  benzyl  phthalate 


1.19 

6.00 

1.32 

7.91 

1.49 

9.23 

1.62 

10.84 

1.56 

9.25 

1.55 

10.02 

1.10 

8.02 

1.30 

9.54 

0.82 

4.11 

1.09 

6.29 

1.59 

10.34 

1.87 

14.02 

3.15 

22.80 

1.14 

6.85 

1.15 

7.68 

1.63 

10.87 

1.12 

6.78 

1.53 

10.06 

2.51 

18.41 

2.23 

16.74 

2.58 

18.82 
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Column:  Zorbax  C-8 
Mobile  phase:  MeOH/water 
Reference:  This  work 

Solute 

Methyl  Paraben 

Ethyl  Paraben 

Propyl  Paraben 

Butyl  Paraben 

Cortisone 

Corticosterone 

Toluene 

Hexylf luorobenzene 

Benzyl  Alcohol 

o-Nitrophenol 

Propachlor 

1-Methyl  naphthalene 

Tri-p-toyly  phosphate 

Methyl  benzyl  amine 

Fluorobenzeze 

Chloropropham 

Dimethyl  phthalate 

Diethyl  phthalate 

Dibutyl  phthalate 

1-Methyl  phenanthrene 

Butyl  benzyl  phthalate 

p-Terphenyl 

Chrysene 


1.11 

5.95 

1.19 

7.10 

1.51 

9.50 

1.66 

10.96 

1.51 

8.60 

1.72 

10.95 

1.23 

8.53 

1.31 

9.07 

0.93 

4.53 

1.07 

6.13 

1.48 

9.40 

1.57 

11.18 

3.41 

24.21 

1.14 

6.65 

1.13 

7.01 

1.65 

11.43 

1.39 

7.92 

1.61 

10.41 

2.73 

19.29 

2.35 

17.31 

2.80 

19.71 

2.82 

20.97 

2.57 

18.19 
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Column:  Zorbax  phenyl 
Mobile  phase:  MeOH/water 
Reference:  This  work 

Solute  Slope  from  equation  3-1     ln-k 

Methyl  Paraben 

Ethyl  Paraben 

Propyl  Paraben 

Butyl  Paraben 

Cortisone 

Corticosterone 

Toluene 

Hexylf luorobenzene 

Benzyl  Alcohol 

o-Nitrophenol 

Propachlor 

1-Methyl  naphthalene 

Methyl  benzyl  amine 

Fluorobenzeze 

Chloropropham 

Dimethyl  phthalate 

Diethyl  phthalate 

Dibutyl  phthalate 

1-Methyl  phenanthrene 

Butyl  benzyl  phthalate 


1.03 

5.50 

1.18 

6.87 

1.49 

8.73 

1.61 

10.41 

1.49 

9.13 

1.69 

11.16 

1.09 

6.79 

1.37 

8.96 

0.80 

3.67 

1.06 

6.00 

1.37 

8.66 

1.81 

12.23 

0.91 

5.20 

0.88 

5.33 

1.60 

10.48 

1.25 

7.39 

1.44 

9.08 

2.62 

18.35 

2.31 

15.78 

2.72 

18.66 
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Column:  Zorbax  TMS 
Mobile  phase:  MeOH/water 
Reference:  This  work 

Solute 

Methyl  Paraben 

Ethyl  Paraben 

Propyl  Paraben 

Butyl  Paraben 

Cortisone 

Corticosterone 

Toluene 

Hexylf luorobenzene 

Benzyl  Alcohol 

o-Nitrophenol 

Propachlor 

1-Methyl  naphthalene 

Tri-p-toyly  phosphate 

Methyl  benzyl  amine 

Fluorobenzeze 

Chloropropham 

Dimethyl  phthalate 

Diethyl  phthalate 

Dibutyl  phthalate 

1-Methyl  phenanthrene 

Butyl  benzyl  phthalate 


In  Kv 


1.11 

5.25 

1.16 

6.11 

1.25 

6.99 

1.56 

8.95 

1.47 

6.97 

1.62 

8.44 

0.91 

4.97 

0.97 

5.61 

0.86 

2.96 

0.93 

4.09 

1.25 

7.09 

1.36 

8.33 

3.18 

20.69 

1.01 

4.94 

0.85 

4.34 

1.44 

8.38 

1.21 

6.26 

1.37 

8.16 

2.21 

14.47 

2.08 

13.27 

2.29 

14.52 
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Column:  Zorbax  cyano 
Mobile  phase:  MeOH/water 
Reference:  This  work 

Solute  Slope  from  equation  3-1      ln_k. 

Methyl  Paraben 

Ethyl  Paraben 

Propyl  Paraben 

Butyl  Paraben 

Cortisone 

Corticosterone 

Toluene 

Hexylf luorobenzene 

Benzyl  Alcohol 

o-Nitrophenol 

Propachlor 

1-Methyl  naphthalene 

Tri-p-toyly  phosphate 

Methyl  benzyl  amine 

Fluorobenzeze 

Chloropropham 

Dimethyl  phthalate 

Diethyl  phthalate 

Dibutyl  phthalate 

1-Methyl  phenanthrene 

Butyl  benzyl  phthalate 


0.87 

3.79 

1.10 

5.12 

1.22 

6.49 

1.35 

7.70 

1.29 

5.81 

1.35 

6.90 

0.86 

4.56 

0.80 

4.52 

0.70 

2.17 

0.74 

3.46 

1.25 

5.99 

1.17 

7.34 

2.50 

15.78 

0.88 

3.83 

0.83 

4.11 

1.26 

7.12 

1.09 

4.99 

1.16 

6.21 

1.84 

11.12 

1.65 

10.83 

1.91 

12.54 
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Column:  Ultrasphere  C-8 
Mobile  phase:  MeOH/water 
Reference:  This  work 

Solute  Slope  from  equation  3-1      In  k- 

Methyl  Paraben 

Ethyl  Paraben 

Propyl  Paraben 

Butyl  Paraben 

Cortisone 

Corticosterone 

Toluene 

Hexylf luorobenzene 

Benzyl  Alcohol 

o-Nitrophenol 

Propachlor 

1-Methyl  naphthalene 

Tri-p-toyly  phosphate 

Methyl  benzyl  amine 

Fluorobenzeze 

Chloropropham 

Dimethyl  phthalate 

Diethyl  phthalate 

1-Methyl  phenanthrene 


1.07 

6.03 

1.27 

7.55 

1.43 

9.35 

1.56 

10.45 

1.44 

8.38 

1.61 

10.31 

1.42 

9.42 

1.53 

10.26 

0.78 

4.03 

0.96 

5.65 

1.40 

8.95 

1.84 

13.12 

3.30 

24.26 

1.02 

5.57 

1.05 

6.97 

1.57 

10.98 

1.25 

7.39 

1.56 

9.85 

2.20 

15.84 

Column:  B  S  J  C-8 
Mobile  phase:  MeOH/water 
Reference:  This  work 
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Solute 

Methyl  Paraben 

Ethyl  Paraben 

Propyl  Paraben 

Butyl  Paraben 

Cortisone 

Corticosterone 

Toluene 

Hexylf luorobenzene 

Benzyl  Alcohol 

o-Nitrophenol 

Propachlor 

1-Methyl  naphthalene 

Tri-p-toyly  phosphate 

Methyl  benzyl  amine 

Fluorobenzeze 

Chloropropham 

Dimethyl  phthalate 

Diethyl  phthalate 

1-Methyl  phenanthrene 


Slope  from  equation  3-1 


In  k. 


1.01 

5.74 

1.12 

6.58 

1.33 

8.44 

1.48 

9.50 

1.37 

7.66 

1.45 

8.91 

0.97 

6.13 

0.99 

6.43 

0.71 

3.24 

0.74 

4.20 

1.17 

7.04 

1.48 

9.80 

3.02 

20.59 

0.86 

4.92 

0.83 

4.96 

1.47 

9.77 

1.10 

6.22 

1.30 

7.84 

1.85 

12.75 
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Column:  B  &  J  C-18 
Mobile  phase:  MeOH/water 
Reference:  This  work 

Solute 

Methyl  Paraben 

Ethyl  Paraben 

Propyl  Paraben 

Butyl  Paraben 

Cortisone 

Corticosterone 

Toluene 

Hexylf luorobenzene 

Benzyl  Alcohol 

o-Nitrophenol 

Propachlor 

1-Methyl  naphthalene 

Methyl  benzyl  amine 

Fluorobenzeze 

Chloropropham 

Dimethyl  phthalate 

Diethyl  phthalate 

1-Methyl  phenanthrene 


ln_k.„ 


1.03 

5.80 

1.17 

7.44 

1.45 

9.11 

1.63 

11.17 

1.53 

9.46 

1.59 

10.34 

1.17 

8.74 

1.26 

9.12 

0.76 

3.97 

0.95 

6.04 

1.40 

9.15 

1.62 

12.14 

1.05 

6.27 

1.12 

7.48 

1.78 

12.34 

1.20 

7.23 

1.46 

9.64 

2.00 

15.22 

Column:  High  density  (Sentell  1987) 
Mobile  phase:  MeOH/water 
Reference:  This  work 
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Solute 

Methyl  Paraben 
Ethyl  Paraben 
Propyl  Paraben 
Butyl  Paraben 
Cortisone 
Corticosterone 
Benzyl  Alcohol 
o-Nitrophenol 
Propachlor 

1-Methyl  naphthalene 
Methyl  benzyl  amine 
Chloropropham 
Dimethyl  phthalate 
Diethyl  phthalate 


3=1 


In-Xv 


0.94 

5.22 

1.05 

6.67 

1.29 

8.41 

1.46 

9.94 

1.34 

6.97 

1.50 

9.87 

0.74 

3.90 

0.75 

4.77 

1.29 

8.41 

1.52 

11.53 

0.82 

4.94 

1.58 

10.91 

1.01 

5.64 

1.36 

9.03 

REFERENCES 

Abel  E.W.,  Pollard  F.H.,  Uden  P.C.,  Nickless  G.;  J.    Chromatogr . ,    1966, 
22,    23. 

Anderson  R.L.;  "Practical  Statistics  for  Analytical  Chemists,"  Van 
Nostrand  Reinhold  Company,  New  York,  1987. 

Antle  P.E.,  Snyder  L.R.;  LC-GC,     1985,  2,    840. 

Antle  P.E.,  Goldberg  A. P.,  Snyder  L.R.;  J.    Chromatogr.,    1985,  321,    1. 

Arai  Y.,  Hirukawa  M.,  Hanai  T.;  J.    Chromatogr.,    1987,  384,    279. 

Baty  J.D.,  Sharp  S.;  J.    Chromatogr.,    1988,  437,  13. 

Bayer  E.,  Paulus  A.;  J.  Chromatogr.,    1987,  400,    1. 

Berendsen  G.E.,  Pikaart  K.A.,  de  Galan  L . ;  J.    Liq.    Chromatogr.,     1980,  3, 
1437. 

Bondi  A.;  J.    Phys .    Chem. ,     1964,  68,     441. 

Braumann  T.;  J.    Chromatogr.,    1986,  313,    191. 

Braumann  T.,  Genieser  H.,  Lullman  C,  Jastorff  B.;  Chromatogr aphia, 
1987,  24,    111. 

Bristow  P. A.;  J.  Chromatogr.,    1978,  149,     13. 

Brown  M.F.;  J.    Chem.    Phys.,    1984,  80,    2808. 

Buszewski  B.,  Suprynowicz  Z.;  Chromatogr  aphia,    1987,  24,    573. 

Buszewski  B.,  Nondek  L.,  Jurasek  A.,  Berek  D . ;  Chromatogr aphia,  1987, 
23,    442. 

Butte  W.,  Fooken  C,  Klussmann  R.,  Schuller  D.;  J.  Chromatogr.,  1981, 
214,     59. 

Cabane  B.;  J.  Phys.    Chem.,    1977,  81,    1639. 

Carr  J.W.,  Harris  J.M.;  Anal.    Chem.,     1986,  58,     626. 

Carr  J.W.,  Harris  J.M.;  A/ial.  Chem.,     1987,  59,     2546. 


123 


124 


Chretien  J.R.,  Walczak  B.,  Allory  L.M.,  Dreux  M.,  Lafosse  M.;  J. 
Chromatogr.,    1986,  371,    253. 

Claudy  P.,  Letoffe  J.M.,  Gaget  C,  Morel  D.,  Serpinet  J.;  J. 
Chromatogr. ,    1985,  329,    331. 

Colin  H.,  Guiochon  G.,  Jandera  P.;  Anal.  Chem. ,     1983,  55,     442. 

Cooper  W.T.,  Lin  L.Y.;  Chromatogr aphia,     1986,  21,    335. 

Delaney  M.F.,  Papas  A.N.,  Walters  M.J.;  J.    Chromatogr. ,    1987,  410,    31. 

Dewaele  C,  Verzele  M.;  J".  Chromatogr.,    1983,  260,     13. 

Dewaele  C,  Mussche  P.,  Verzele  M.;  HRC&CC,    1982,  5,     616. 

Dill  K.A.;  J.    Phys.    Chem.,    1987,  91,    1980. 

Dill  K.A.,  Flory  P . J . ;  Proc.    Natl.    Acad.    Sci.,     1980,  77,  3115. 

Dill  K.A.,  Flory  P . J . ;  Proc.    Natl.    Acad.    Sci.,     1981,  78,     676. 

Dill  K.  A.,  Koppel  D.E.,  Cantor  R.S.,  Dill  J.D.,  Bendedouch  D.,  Chen 
S.H.;  ;  Nature,  1984,  309,    42. 

Dill  K.A.,  Naghizadeh  J.,  Marqusee  J. A.;  Ann.  Rev.    Phys.    Chem.,    1988, 
39,    425. 

Dolan  J.W.,  Gant  J.R.,  Snyder  L.R.;  J.    Chromatogr. ,    1979,  165,    31. 

Dorsey  J.G.,  Dill  K.A.;  Chem.    Rev.,     1989,  89,     331. 

Dorsey  J.G.,  Johnson  B.P.;  J.    Liq.    Chromatogr. ,    1987,  10,    2695. 

Engelhardt  H.,  Dreyer  B.,  Schmidt  H . ;  Chromatogr aphia,    1982,  16,    11. 

Engelhardt  H.,  Miiller  H.;  J.    Chromatogr.,    1981,  218,    395. 

Engelhardt  H.,  Miiller  H.,  Dreyer  B.;  Chromatographia,     1984,  19,    240. 

Evans  M.B.,  Dale  A.D.,  Little  C.J.;  Chromatographia,     1980,  13,     5. 

Feng  Y.,  Zhu  P.,  Hu  Z  .  ;  Chromatographia,     1988,  25,     382. 

Funasaki  N.,  Hada  S.,  Neya  S.;  J.    Chromatogr. ,    1986,  361,    33. 

Gazda  K.,  Kaminski  M.,  Klawiter  J.,  Kowalczyk  J.S.,  Makuch  B., 
Prusiewicz  K.,  Sledzinska  B.;  J.    Chromatogr. ,    1980,  191,     9. 

Gilpin  R.K.;  J.    Chromatogr.    Sci.,    1984,  22,    371. 

Gilpin  R.K.,  Gangoda  M.E.;  J.    Chromatogr.    Sci.,    1983,  21,    352. 

Gilpin  R.K.,  Gangoda  M.E.;  Anal.  Chem.,    1984,  56,     1470. 


125 


Gobet  J.,  Kovats  E.;  Adsorption    Sci.    &    Technol.,    1984,  1,  77. 

Goldberg  A. P.;  Anal.    Chem.,    1982,  54,    342. 

Golding  R.D.,  Barry  A.J.,  Burke  M.F.;  J.    Chromatogr. ,     1987,  384,     105. 

Hammers  W.E.,  Meurs  G.J.,  Lingy  C.L.;  J.    Chromatogr. ,     1982,  247,  1. 

Hanai  T.,  Hubert  J.;  HRC&CC,     1983,  6,    20. 

Hanai  T . ,  Hubert  J.;  J.    Chromatogr.,    1984a,  290,    197. 

Hanai  T.,  Hubert  J.;  J.    Chromatogr.,    1984b,  302,    89. 

Hansen  S.H.,  Helboe  P.,  Thomsen  M.;  J.    Chromatogr. ,  1986,  363,    39. 

Hansen  S.H.,  Helboe  P.,  Thomsen  M.;  Tr.    Anal.    Chem.,    1988,  7,  389. 

Hennion  M.C.,  Rosset  R. ;  Chromatogaphia . ,  1988,  25,    43. 

Hildebrand  J.H.,  Scott  R.L.;  "The  Solubility  of  Nonelectrolytes, " 
Reinhold,  New  York,  1950. 

Hill  T.L.;  "An  Introduction  to  Statistical  Thermodynamics,"  Dover 
Publications  Inc.,  New  York,  1986. 

Horvath  C.,  Melander  W.,  Molnar  I.;  J.    Chromatogr. ,  1976,  125,     129. 

Howard  G.A.,  Martin  A. J. P.;  Biochem.    J.,     1950,  46,     532. 

Jaroniec  M.;  J.    Liq.    Chromatogr. ,    1984,  suppl.2,    393. 

Jandera  P.;  Chromatographia,    1985,  19,    101. 

Jandera  P.,  Colin  H.,  Guiochon  G.;  Anal.  Chem.,    1982,  54,     435. 

Jinno  K.,  Kawasaki  K.;  Chromatographia,    1983a,  17,  337. 

Jinno  K.,  Kawasaki  K . ;  Chromatographia,     1983b,  17,  445. 

Johnson  B.P.,  Khaledi  M.G.,  Dorsey  J.G.;  Anal.  Chem.;    1986,  58,     2354. 

Johnson  B.P.;  "Solvatochromic  Solvent  Polarity  Measurements,  Retention, 
and  Selectivity  in  Reversed  Phase  Liquid  Chromatography",  Doctoral 
Dissertation,  University  of  Florida,  1986. 

Kaliszan  R.;  CRC  Crit .    Rev.    Anal.    Chem.,     1986,  16,     323. 

Kamlet  M.J.,  Abboud  J.M.,  Abraham  M.H.,  Taft  R.W.;  J.  Org.    Chem.,     1983, 
48,    2877. 


126 

Karger  B.L.,  LePage  J.N.,  Tanaka  N . ;  "High  Performance  Liquid 
Chromatography- Advances  and  Propectives,  Vol.  2,"  Academic  Press,  New 
York,  1980,  p. 133. 

Karger  3.L.,  Snyder  L.R.,  Eon  C;  Anal.  Chem.  ,  1978,  50,     2126. 

Khong  T.M.,  Simpson  C . F . ;  Chromatographia,    1987,  24,    385. 

Kirkland  J.J.,  DeStefano  J . J . ;  J.    Chromatogr.    Sci.,    1970,  8,    309. 

Kosower  E.M.;  J.  Am.  Chem.    Soc,     1958,  80,     3253. 

Kovats  E.;  "Advances  in  Chromatography,  Vol.  1,"  Marcel  Dekker,  New 
York,  1965. 

Lipford  L.C.;  "Temperature  and  Mobile  Phase  Composition  Effects  on  the 
Retention  Behavior  of  Polynuclear  Aromatic  Hydrocarbons:  Thermodynamics 
and  Stationary  Phase  Formation  Considerations,"  Masters  Thesis, 
University  of  Florida,  1985. 

Lochmiiller  C.H.,  Hangac  H.H.,  Wilder  D.R.;  J.    Chromatogr.    Sci.,    1981, 
19,    130. 

Lochmiiller  C.H.,  Hunnicutt  M.L.,  Mullaney  J .  F .  ;  J.    Phys .    Chem.,     1985, 
89,    5770. 

Lochmiiller  C.H.,  Marshall  D.B.;  Anal.  Chim.    Acta;    1982,  142,     63. 

Lochmiiller  C.H.,  Wilder  D.R.;  J.    Chromatogr.    Sci.,    1979,  17,  574. 

Lork  K.D.,  Unger  K.K.,  Kinkel  J.N.;  J.    Chromatogr.,     1986,  352,     199. 

Majors  R.E.;  "High  Performance  Liquid  Chromatography-Advances  and 
Propectives,  Vol.  1,"  Academic  Press,  New  York,  1980,  p. 78. 

Marqusee  J. A.,  Dill  K.A.;  J.    Chem.    Phys.,     1986a,  85,     434. 

Marqusee  J. A.,  Dill  K.A.;  Macromolecules,     1986b,  19,    2420. 

Marshall  D.3.,  Cole  C.L.,  Connolly  D.E.;  J.    Chromatogr.,    1986,  361,    71. 

Marshall  D.3.,  Cole  C.L.,  Norman  A.D.;  J.    Chromatogr.    Sci.,    1987,  25, 
262. 

Marshall  D.3.,  Stutler  K.A.,  Lochmiiller  C.H.;  J.    Chromatogr.    Sci.,     1984, 
22,    217. 

Martire  D.E.,  Boehm  R.E.;  J.    Phys.    Chem.,     1983,  81,     1045. 

McCormick  R.M.,  Karger  B.L.;  Anal.  Chem.,    1980a,  52,    2249. 

McCormick  R.M.,  Karger  B.L.;  J.    Chromatogr. ,     1980b,  199,     259. 

McNally  M.E.,  Rogers  L.B.,;  J.    Chromatogr. ,     1985,  331,    23. 


127 


McReynolds  W.O.;  J.    Chromatogr.    Sci.,    1970,  8,     685. 

Melander  W.,  Horvath  C . ;  "High  Performance  Liquid  Chromatography- 
Advances  and  Propectives,  Vol.  1",  Academic  Press,  New  York,  1980, 
p. 113. 

Melander  W.R.,  Erard  J.F.,  Horvath  C .  ;  J.    Chromatogr.,    1983,  282,    211. 

Melander  W.R.,  Stoveken  J.,  Horvath  C . ;  J.    Chromatogr. ,  1980,  199,    35. 

Michels  J. J.;  "Solvatochromic  Investigation  of  Chromatographic 
Processes,"  Doctoral  Dissertation,  University  of  Florida,  1989. 

Minick  D.J.,  Sabatka  J.J.,  Brent  D.A.,;  J.    Liq.    Chromatogr. ,  1987,  10, 
2565. 

Miyake  K.,  Mizuno  N.,  Terada  H.;  J.    Chromatogr. ,     1988,  439,    227. 

Mockel  H.J.,  Welter  G.,  Meltzer  H.;  J.    Chromatogr. ,  1987,  388,    255. 

Nawrocki  J.;  J.    Chromatogr. ,    1986,  362,     117. 

Nawrocki  J.;  J.    Chromatogr. ,    1987,  391,    266. 

Nawrocki  J.,  Buszewski  B.;  J.    Chromatogr. ,    1988,  449,     1. 

Onsager  L . ;  J.    Am.    Chem.    Soc,    1936,  58,    1486. 

Opperhuizen  A.,  Sinnige  T.L.,  van  der  Steen  J.M.D.,  Hutzinger  O.;  J. 
Chromatogr.,    1987,  388,    51. 

Reymond  D.,  Chung  G.N.,  Mayer  J.M.,  Testa  B.;  J.    Chromatogr. ,    1987,  391, 
97. 

Rohrschneider  L.;  J.    Chromatogr. ,    1965,  17,    1. 

Sadek  P.C.,  Carr  P.W.;  J".  Chromatogr.    Sci.,     1983,  21,    314. 

Sadek  P.C.,  Carr  P.W.;  J.    Chromatogr.,    1984,  288,    25. 

Sadek  P.C.,  Carr  P.W.,  Doherty  R.M.,  Kamlet  M.J.,  Taft  R.W.,  Abraham 
M.H.;  Anal.  Chem.,    1985,  57,  2971. 

Sadek  P.C.,  Koester  C.J.,  Bowers  L.D.;  J.    Chroamtogr.    Sci.,     1987,  25, 
489. 

Sander  L.C.;  J.  Chroamtogr.    Sci.,    1988,  26,    380. 

Sander  L.C.,  Callis  J.B.,  Field  L.R.;  Anal.  Chem.,    1983,  55,     1068. 

Sander  L.C.,  Field  L.R.;  Anal.  Chem.,    1980,  52,    2009. 

Sander  L.C.,  Wise  S.A.;  Anal.  Chem.,    1984,  56,    504. 


128 


Sander  L.C.,  Wise  S.A.;  HRC&CC,    1988,  11,    383. 

Schantz  M.,  Martire  D.E.;  J.    Chromatogr . ,     1987,  391,    35. 

Schoenmakers  P.J.,  Billiet  H.A.H.,  de  Galan  L.J.;  J.    Chromatogr. ,    1979, 
185,     179. 

Schoenmakers  P.J.,  Billiet  H.A.H.,  de  Galan  L.J.;  J.    Chromatogr.,    1981, 
205,    13. 

Schoenmakers  P.J.,  Billiet  H.A.H.,  de  Galan  L.J.;  Chromatographia,    1982, 
15,    205. 

Schoenmakers  P.J.,  Billiet  H.A.H.,  de  Galan  L.J.;  J.    Chromatogr.,    1983, 
282,     107. 

Schoenmakers  P.J.,  Billiet  H.A.H.,  Tijssen  R.,  de  Galan  L.J.;  J. 
Chromatogr.,    1978,  149,    519. 

Sentell  K.B.;  "Interphase  Solubility  and  Chromatographic  Retention," 
Doctoral  Dissertation,  University  of  Florida,  1987. 

Sentell  K.B.,  Barnes  K.W.,  Dorsey  J.G.;  J.    Chromatogr.,    1988,  455,     95. 

Sentell  K.B.,  Dorsey  J.G.;  Anal.  Chem. ,    1989a,  61,     930. 

Sentell  K.B.,  Dorsey  J.G.;  J.    Chromatogr.,    1989b,  461,    193. 

Sharaf  M.A.,  Illman  D.L.,  Kowalski  B.R.;  "Chemometric" ,  John  Wiley  & 
Sons,  New  York,  1986. 

Sinanoglu  O . ;  Chem.    Phys .    Lett.,    1967,  1,    340. 

Slaats  E.H.,  Markovski  W.,  Fekete  J.,  Poppe  H.;  J.    Chromatogr. ,     1981, 
207,  299. 

Smith  R.J.,  Nieass  C.S.,  Wainwright  W.S.;  J.    Liq.    Chromatogr. ,    1986,  9, 
1387. 

Smith  R.M.;  J.    Chromatogr.,    1981,  209,    1. 

Smith  R.M.;  J.    Chromatogr. ,    1982a,  236,    313. 

Smith  R.M.;  J.    Chromatogr. ,     1982b,  236,     321. 

Smith  R.M.,  Hurdley  T.G.,  Gill  R.,  Moffat  A.C.;  J.    Chromatogr. ,     1986, 
351,    259. 

Smith  R.M.,  Miller  S . L . ;  J.    Chromatogr.,     1989,  464,    297. 
Snyder  L.R.;  J.    Chromatogr. ,    1974,  92,    233. 
Snyder  L.R.;  J.    Chromatogr.    Sci.,    1978,  16,    233. 


129 


Snyder  L.R.,  Dolan  J.W.,  Gant  J.R.;  J.    Chromatogr. ,    1979,  165,    3. 

Stewart  H.N.M.,  Perry  S.G.;  J.    Chromatogr. ,     1968,  37,     97. 

Tanaka  N.,  Sakagami  K.,  Araki  M.;  J.  Chromatogr. ,    1980,  199,    327. 

Tanaka  N.,  Tokuda  Y.,  Iwaguchi  K.f  Araki  M . ;  J.    Chromatogr. ,     1982,  239, 
761. 

Tanford  C . ;  "The  Hydrophobic  Effect,  2nd  ed.,"  Wiley-Intersciences,  New 
York,  1980. 

Taft  R.W.,  Abboud  J.M.,  Kamlet  M.J.,  Abraham  M.H.;  J.    Solution   Chem., 
1985,  14,    153. 

Taft  R.W.,  Kamlet  M.J.;  J".  Am.  Chem.  Soc,     1976,  98,     2886. 

Tijssen  R.,  Billiet  H.A.H.,  Schoenmakers  P.J.;  J.    Chromatogr. ,    1976, 
122,    185. 

Unger  K.  K.;  "Porous  Silica,"  Elservier  Scientific  Publishing  Company, 
New  York,  1979. 

Verzele  M.,  LC-Mag.    Liq.    Chromatogr.    HPLC,     1983,  J,  217. 

Verzele  M.,  DePotter  M.,  Ghysels  J.;  HRC&CC,     1979,  2,  151. 

Verzele  M.,  Dewaele  C;  J.  Chromatogr. ,    1981,  217,  399. 

Verzele  M.,  Dewaele  C.,  Duquet  D .  ;  J".  Chromatogr. ,    1985,  329,    351. 

Wahlund  K.G.,  Sokolowski  A.;  J.    Chromatogr. ,    1978,  151,    299. 

Walczak  B.,  Allory  L.M.,  Lafosse  M.,  Dreux  M.,  Chretien  J.R.;  J. 
Chromatogr.,    1987,  395,    183. 

Wallings  C;  J.    Am.    Chem.    Soc,    1970,  72,  1164. 

Wells  M.J.M.,  Clark  C.R.,  Patterson  R.M.;  J.    Chromatogr. ,    1982,  235,    61. 

Wells  M.J.M.,  Clark  C.R.,  Patterson  R.M.;  Anal.  Chem.,    1986,  58,    1625. 

White  S.H.,  King  G.I.,  Cain  J.E.;  Nature,     1981,  290,     161. 

Wise  S.A.,  Sander  L.C.;  HRC&CC,    1985,  8,    248. 

Woodburn  K.B.;  "Thermodynamics  and  Mechanisms  of  Sorption  for 
Hydrophobic  Organic  Compounds  on  Natural  and  Artifical  Sorbent 
Materials,"  Doctoral  Dissertation,  University  of  Florida,  1985. 

Yonker  C.R.,  Zwier  T.A.,  Burke  M.F.;  J.    Chromatogr. ,  1982a,  241,    257. 

Yonker  C.R.,  Zwier  T.A.,  Burke  M . F . ;  J.    Chromatogr . ,     1982b,  241,    269. 


BIOGRAPHICAL  SKETCH 
Peter  Tai  Yuen  Ying  was  born  on  February  28,  1962,  in  Hong  Kong. 
He  completed  his  high  school  education  in  Hong  Kong,  and  at  the  age  of 
19,  he  came  over  to  America  to  further  his  education.   He  first  attended 
the  University  of  Wisconsin-Madison.   Due  to  numerous  reasons,  and  one 
being  the  dislike  of  long  winters,  he  transferred  to  the  College  of 
Charleston  in  South  Carolina  in  June  1983.   He  later  graduated  with  a 
B.S.  degree  in  chemistry  in  May  1985  from  the  College  of  Charleston.   He 
entered  the  graduate  school  at  the  University  of  Florida  in  August  1985. 
After  all  these  years  of  living  and  learning  in  America,  he  has  become 
fascinated  by  the  western  cuisines,  music  and  women  (but  not  necessarily 
in  that  order) .   Therefore,  upon  completion  of  the  requirements  for  the 
degree  of  Doctor  of  Philosophy,  he  will  join  the  Wesley-Jessen  Corp.  in 
Chicago  as  a  research  chemist  so  he  can  continue  to  pursue  his 
fascinations . 


130 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it 
conforms  to  acceptable  standards  of  scholarly  presentation  and  is  fully 
adequate,  in  scope  and  quality,  as  a  dissertation  for  the  degree  of 
Doctor  of  Philosophy. 


John  G\  Dorsey,  Chair 

Associate  Professor  of  Chemistry 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it 
conforms  to  acceptable  standards  of  scholarly  presentation  and  is  fully 
adequate,  in  scope  and  quality,  as  a  dissertation  for  the  degree  of 
Doctor  of  Philosophy. 


h^^UUjir^/xJyy^ 


/ames  D.  Winefordner/ 

Graduate  Research  p/ofessor  of  Chemistry 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it 
conforms  to  acceptable  standards  of  scholarly  presentation  and  is  fully 
adequate,  in  scope  and  quality,  as  a  dissertation  for  the  degree  of 
Doctor  of  Philosophy. 


Associate  Professor  of  Chemistry 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it 
conforms  to  acceptable  standards  of  scholarly  presentation  and  is  fully 
adequate,  in  scope  and  quality,  as  a  dissertation  for  the  degree  of 
Doctor  of  Philosophy. 


u 


k'V\£.> 


/   N.MfcM, 


T-<  1 


Vaneica  Y.  Young  f!  ^ 

Associate  Professor  of  Chejnistry 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it 
conforms  to  acceptable  standards  of  scholarly  presentation  and  is  fully 
adequate,  in  scope  and  quality,  as  a  dissertation  for  the  degree  of 
Doctor  of  Philosophy. 


;s4e  F.  Gregory      (^J     ^J 


Jesrse  F.  Grec 

Professor  of  Food  Science  &  Human 

Nutrition 


This  dissertation  was  submitted  to  the  Graduate  Faculty  of  the 
Department  of  Chemistry  in  the  College  of  Liberal  Arts  and  Sciences  and 
to  the  Graduate  School  and  was  accepted  as  partial  fulfillment  of  the 
requiremnets  for  the  degree  of  Doctor  of  Philosophy. 


August,  1989 


Dean,  Graduate  School 


UNIVERSITY  OF  FLORIDA 

3  1262  08553  6414 


